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We study large-distance contributions to the nucleon's parton densities in the transverse coordi- 
nate (impact parameter) representation based on generalized parton distributions (GPDs). Chiral 
dynamics generates a distinct component of the partonic structure, located at momentum fractions 
X < ALtt/Mn and transverse distances h ~ l/M-n- We calculate this component using phenomenolog- 
ical pion exchange with a physical lower limit in b (the transverse "core" radius estimated from the 
nucleon's axial form factor, -Rcoro = 0.55 fm) and demonstrate its universal character. This formula- 
tion preserves the basic picture of the "pion cloud" model of the nucleon's sea quark distributions, 
while restricting its application to the region actually governed by chiral dynamics. It is found that 
(a) the large-distance component accounts for only ~ 1/3 of the measured antiquark flavor asym- 
metry d — M at a; ~ 0.1; (b) the strange sea quarks, s and s, are significantly more localized than 
the light antiquark sea; (c) the nucleon's singlet quark size for x < 0.1 is larger than its gluonic size, 
q+q > ^ suggested by the t-slopes of deeply-virtual Compton scattering and exclusive 

J/ if) production measured at HERA and FNAL. We show that our approach reproduces the general 
A^c-scaling of parton densities in QCD, thanks to the degeneracy of A'^ and A intermediate states 
in the large- A^c limit. We also comment on the role of pionic configurations at large longitudinal 
distances and the limits of their applicability at small x. 
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I. INTRODUCTION 

Parton densities summarize the structure of the nu- 
cleon probed in high-momentum transfer processes such 
as deep-inelastic lepton-nucleon scattering and produc- 
tion of high-mass systems (jets, heavy particles) in 
nucleon-nucleon collisions. They are defined in the con- 
text of a factorization procedure, by which the cross sec- 
tion of these processes is separated into a short-distance 
quark/gluon subprocess, calculable in perturbative QCD, 
and the distribution of the partons in the initial state, 
and thus represent long-distance, low-energy character- 
istics of the nucleon. As such, they are governed by 
the same low-energy dynamics which determines other 
nucleon observables like the vector and axial couplings 
(to which they are related by the partonic sum rules), 
form factors, meson-nucleon couplings etc. Of partic- 
ular interest are the charge [u — u,d — d) and flavor 
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(u — d,u — d) non-singlet quark densities, which exhibit 
only weak scale dependence and are of non-perturbative 
origin; they represent quasi-observables which directly 
probes the QCD quark structure of the nucleon at low 
resolution scales. 

The long-distance behavior of strong interactions at 
low energies is governed by the spontaneous breaking of 
chiral symmetry in QCD. The Goldstone boson nature 
of the pion explains its small mass on the hadronic scale 
and requires its coupling to other hadrons to vanish in 
the long- wavelength limit. The resulting "chiral dynam- 
ics" gives rise to a number of distinctive phenomena at 
distance scales ^ such as the tttt, ttA^ and NN 

interaction at large distances, the pion pole in the axial 
current matrix element, etc. An important question is 
how chiral dynamics affects the nucleon's parton densi- 
ties, and whether one can see any signs of chiral effects 
in observables of high-momentum transfer processes. 

The prime candidate for an effect of chiral dynam- 
ics in parton densities has been the flavor asymmetry 
of the light antiquark densities in the nucleon. Mea- 
surements of the proton-neutron structure function dif- 
ference in inclusive deep-inelastic scattering [l|, semi- 
inclusive meson production and particularly Drell- 
Yan pair production 0, Q have unambiguously shown 
that [d — u] (x) > in the proton for x < 0.3, and have 
partly mapped the x-dependence of the asymmetry; see 
Refs. Ji for a review of earlier experimental results. The 
basic picture is that the "bare" proton can make a transi- 
tion to a virtual state containing a pion, and fluctuations 
p n7r+ are more likely than p — > A++7r^, resulting in 
an excess of tt"*" over 7r~ in the "dressed" proton. Fol- 
lowing the original prediction of Ref . [Tj] , which included 
only the nucleon intermediate state, this idea was imple- 
mented in a variety of dynamical models, which incorpo- 
rate finite-size effects through various types of hadronic 
form factors associated with the ttNN and ttNA ver- 
tices; see Refs. @ for a review of the extensive litera- 
ture. It was noted long ago Q that in order to reproduce 
the fast decrease of the observed asymmetry with x one 
needs ttNN form factors much softer than those com- 
monly used in meson exchange parametrizations of the 
NN interaction However, even with such soft form 
factors the pion transverse momenta in the nucleon gen- 
erally extend up to values ^ This raises the 
question to what extent such models actually describe 
long-distance effects associated with soft pion exchange 
(momenta ~ M^), and what part of their predictions 
is simply a parametrization of short-distance dynamics 
which should more naturally be described in terms of 
non-hadronic degrees of freedom. More generally, one 
faces the question how to formulate the concept of the 
"pion cloud" in the nucleon's partonic structure [ll| in a 
manner consistent with chiral dynamics in QCD. 

A framework which allows one to address these ques- 
tions in a systematic fashion is the transverse coordi- 
nate (or impact parameter) representation, in which the 
distribution of partons is studied as a function of the 



longitudinal momentum fraction, x, and the transverse 
distance, 6, of the parton from the transverse center-of- 
momentum of the nucleon p^.[T3|. In this representation, 
chiral dynamics can be associated with a distinct compo- 
nent of the partonic structure, located at x ^ M^/Mat 
and b ^ l/M^r. In a previous work (T^ we have shown 
that in the gluon density this large-distance component 
is sizable and causes the nucleon's average gluonic trans- 
verse size, to grow if x drops below Mt^/M^, in 
agreement with the t-slopes observed in exclusive J/ip 
photo- and electroproduction at HERA [H, ^ , FNAL 
|17| . and experiments at lower energies. Essential in this 
is the fact that in the gluon density (more generally, in 
any isoscalar parton density) the pion cloud contribu- 
tions from N and A intermediate states have the same 
sign and add constructively. The special role of the A 
compared to other excited baryon states is supported by 
the fact that in the large-A^c limit of QCD the N and A 
are degenerate and enter on an equal footing. 

In this article we perform a comprehensive study of the 
chiral large-distance component of the nucleon's partonic 
structure, considering both its contribution to the total 
quark/antiquark/gluon densities and to the nucleon's av- 
erage partonic transverse size. The method we use to 
calculate this component is phenomenological pion ex- 
change formulated in the impact parameter representa- 
tion, restricted to the region of large transverse distances. 
A physical lower limit in b for ttB [B = N, A) configu- 
rations in the nucleon wave function is set by the trans- 
verse "core" radius, estimated from the nucleon's axial 
form factor, i?corc — 0.55 fm, and we explicitly demon- 
strate the universal character of the pionic contributions 
in the region b > -Rcoro- This formulation preserves the 
basic physical picture of the "pion cloud" model of the 
nucleon's sea quark distributions, while restricting its ap- 
plication to the region actually governed by chiral dy- 
namics. In fact, our study serves both a conceptual and 
a practical purpose. First, we want to establish in which 
region of transverse distances the results of the tradi- 
tional pion cloud model are model-independent and can 
be associated with large-distance chiral dynamics. Sec- 
ond, we want to employ this model to actually calcu- 
late the universal large-distance component and study 
its properties. A preliminary account of our study of the 
fiavor asymmetry d — u was presented in Ref. [l8| . 

The investigation reported here proceeds in several 
steps. In Sec. [TTl we develop the theory of large-distance 
contributions to the partonic structure from a general, 
model-independent perspective. We outline the para- 
metric region of ttB configurations in the nucleon wave 
function, the properties of the ^-dependent momentum 
distribution of pions in the nucleon, its large-5 asymp- 
totics, and the convolution formulas for the nucleon par- 
ton densities. In Sec. lIIIi we investigate the phenomeno- 
logical pion cloud model in the impact parameter repre- 
sentation, and demonstrate that at large b its predictions 
become independent of the t:NB form factors modeling 
the short-distance dynamics. We also comment on the 
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extension of this model to SU{3) flavor. In Sec. IIVI we 
then apply this model to calculate the large-distance con- 
tributions to the sea quark distributions in the nucleon, 
including the isovector (d — u) and isoscalar (u + d) light 
quark sea, the strange sea (s,s), and the S'f7(3)-flavor 
symmetry breaking asymmetry (u + d — 2s). We com- 
pare the calculated large-distance contributions to em- 
pirical parametrizations of the parton densities and thus 
indirectly infer the contribution from the short-distance 
region ("core"), which cannot be calculated in a model- 
independent way. In the course of this we see how the 
restriction to large b solves several problems inherent in 
the traditional pion cloud model which formally allows 
for pionic configurations also at small impact parameters. 
In Sec.|V]we consider the large-distance contributions to 
the nucleon's partonic transverse size (6^), which is acces- 
sible experimentally through the ^-slope of hard exclusive 
processes j*N ^ M + N {M = meson, 7, etc.). Because 
of the emphasis on large distances this quantity is cal- 
culable in a practically model-independent manner and 
represents a clean probe of chiral dynamics in the par- 
tonic structure. Specifically, we show that at a: ~ 10~^ 
the large-distance contribution to the nucleon's singlet 
quark transverse size, is larger than that to the 

gluonic size, which is consistent with the observa- 

tion of a larger t-slope in deeply-virtual Compton scat- 
teri ng il9, . ,2Qj than in exclusive J /ip production at HERA 
[li|l|~In SecEI] we discuss the correspondence of the 
phenomenological pion exchange contribution to the nu- 
cleon parton densities with the large- A'c limit of QCD. 
In particular, we show that the large-distance contribu- 
tions obtained from pion exchange reproduce the general 
A^c-scaling of parton densities in QCD, thanks to the de- 
generacy of N and A intermediate states in the large- A'c 
limit. This re-affirms the need to include intermediate 
A states on the same footing as the nucleon, and shows 
that the phenomenological large-distance contributions 
considered here are a legitimate part of the nucleons par- 
tonic structure in large-A^c QCD. Finally, in Sec. IVIII 
we focus on the physical limitations to the picture of in- 
dividual -kB configurations at small x, arising from the 
non-chiral growth of the transverse sizes due to diffu- 
sion, and from chiral corrections to the structure of the 
pion. We also comment on the role of chiral dynamics at 
large longitudinal distances. Our summary and outlook 
are presented in Sec. IVIIII The two appendices present 
technical material related to the meson-nucleon coupling 
constants for SU (3) flavor symmetry, and the numerical 
evaluation of the 6-dcpendcnt pion momentum distribu- 
tions in the nucleon. 

In the context of our studies of the strange sea 
quark distributions, s{x) and s{x), and the SU{2i) fla- 
vor symmetry-breaking asymmetry, \u -\- d — 2s] (x) , we 
consider also contributions from configurations contain- 
ing SU[i) octet mesons {KA,KT,,KT.* .rjN) to the nu- 
cleon's partonic structure at large distances. While such 
high-mass configurations are not governed by chiral dy- 
namics and treated at a purely phenomenological level, it 



is interesting to compare their large-distance tails with 
those of chiral contributions from pions. We note that the 
issue of the strange sea in the nucleon (s, s) and the ques- 
tion of possibly different x-distributions of s and s has 
acquired new urgency following the results of the NuTeV 
experiment in semi-inclusive charged-current neutrino 
DIS, which can discriminate between s and s via the pro- 
cess W+ + s^c^^. 

Chiral contributions to the nucleon's parton densities 
have been studied extensively within chiral perturbation 
theory (23l . [23 | , mostly with the aim of extrapolating lat- 
tice QCD results obtained at large pion masses toward 
lower values [2^. Chiral perturbation theory was also 
applied to GPDs, including the impact parameter rep- 
resentation [2^ [23, [2I, [2§| . Compared to these calcula- 
tions, which use methods of effective field theory based 
on a power-counting scheme, we take here a more prag- 
matic approach. We study the pion distribution in the 
nucleon in a phenomenological approach which incorpo- 
rates the finite bare nucleon size through form factors, 
and investigate numerically in which region the results 
become insensitive to the form factors and can be at- 
tributed to universal chiral dynamics [6^ . In this ap- 
proach we maintain exact relativistic kinematics (physi- 
cal pion and nucleon masses) and calculate distributions 
of finite support, which are then analyzed in the differ- 
ent parametric regions and matched with the asymptotic 
"chiral" predictions. This also allows us to deal with 
the strong cancellations between contributions from N 
and A intermediate states in the isovector quark densi- 
ties, which are difficult to accommodate within a power 
counting scheme. In fact, the cancellation becomes ex- 
act in the large-iVc limit of QCD and ensures the proper 
1 /Nc counting required of the isovector antiquark distri- 
bution in QCD dl]. 

In this study we focus on chiral large-distance contri- 
butions to the nucleon's partonic structure at moderately 
small momentum fractions, x 10~^, which arise from 
individual ttB [B = N, A) configurations in the nucleon 
wave function. When extending the discussion toward 
smaller x, several effects need to be taken into account 
which potentially modify this picture. One is diffusion in 
the partonic wave function, which causes the transverse 
size of the nucleon's partonic configurations to grow at 
small X (however, this effect is suppressed at large Q"^). 
Another effect are possible chiral corrections to the struc- 
ture of the pion itself, which were recently studied in an 
approach based on resummation of chiral perturbation 
theory in the leading logarithmic approximation 31]. We 
discuss the limitations to the applicability of the picture 
of individual ttB configurations in Sees . I VlfXl and fVII Bl 
We also comment on the role of nB configurations at 
large longitudinal separations and arbitrary transverse 
distances, and point out that there may be a window for 
a chiral regime at x 10~^; at smaller x coherence 
effects become dominant; see Sec. I VII CI A detailed in- 
vestigation of this new regime will be the subject of a 
separate study. 



4 



II. CHIRAL DYNAMICS AND PARTONIC 
STRUCTURE 

A. Parametric region of chiral component 

As the first step of our study we want to delineate the 
parametric region where parton densities are governed 
by chiral dynamics and establish its numerical limits, as 
imposed by other, non-chiral physical scales. The pri- 
mary object of our discussion is the pion longitudinal 
momentum and transverse coordinate distribution in a 
fast-moving nucleon, fniy,b), where y is the pion mo- 
mentum fraction. Here we introduce this concept heuris- 
tically, appealing to its obvious physical meaning; its pre- 
cise definition in terms of GPDs and its region of appli- 
cability will be elaborated in the following. 

Chiral dynamics generally governs contributions to nu- 
cleon observables from large distances, of the order 1 /M^, 
which is assumed here to be much larger than all other 
hadronic length scales in question. These contributions 
result from exchange of "soft" pions in the nucleon rest 
frame; in the time-ordered formulation of rclativistic dy- 
namics these are pions with energies E-j^ ^ M-^ and mo- 
menta Ifcjrl ~ -Mtt- Chiral symmetry provides that such 
pions couple weakly to the nucleon and to each other, so 
that their eflFects can be computed perturbatively. Boost- 
ing these weakly interacting pion-nucleon configurations 
to a frame in which the nucleon is moving with large ve- 
locity, we find that they correspond to longitudinal pion 
momentum fractions of the order [63| 

y ^ M^/Mn. (1) 

At the same time, the soft pions' transverse momenta, 
which are not affected by the boost, correspond to trans- 
verse distances of the order 

b ^ l/M^. (2) 

Together, Eqs. ^ and ^ determine the parametric re- 
gion where the pion distribution in the fast-moving nu- 
cleon is governed by chiral dynamics, and the soft pion 
can be regarded as a "parton" in the nucleon's wave func- 
tion in the usual sense (see Fig.[T]). 

The condition Eq. ([1]) implies that the pion momentum 
fraction in the nucleon is parametrically small, y \, 
i.e., the soft pion is a "slow" parton. As a consequence, 
one can generally neglect the recoil of the spectator sys- 
tem and identify the distance b with the separation of 
the pion from the transverse center-of-momentum of the 
spectator system, r = b/{l — y) [g^l [13 • This circum- 
stance greatly simplifies the spatial interpretation of chi- 
ral contributions to the parton densities. 

Pionic configurations in the nucleon wave function arc 
physically meaningful only if the transverse separation 
of the pion and the spectator system is larger than the 
sum of the intrinsic "non-chiral" sizes of these objects. 
This basic fact imposes a limit on the applicability of 
chiral dynamics, even though the dynamics itself may not 




FIG. 1: Parametric region where the pion distribution in the 
nucleon is governed by chiral dynamics. The variables are 
the pion longitudinal momentum fraction, y, and transverse 
position, b. 



change dramatically at the limiting distance. In order 
to make the picture of Fig. [1] quantitative, we have to 
estimate down to which values of b the concept of pionic 
configurations is applicable. 

The transverse size of the "core" in the nucleon's par- 
tonic wave function in the valence region {x ^ 10^^) can 
be estimated from the transverse axial charge radius of 
the nucleon, which does not receive contributions from 
the pion cloud Hi, HI]: 

(6')axial = |(r2)axial ~ 0.3 fm^, (3) 

where the factor 2/3 results from converting the 3- 
dimensional charge radius in the rest frame into the 2- 
dimensional transverse charge radius in the frame where 
the nucleon is moving fast. Identifying the core radius 
with the transverse RMS radius, we obtain 

Rcorc = [ (6') axial] ~ 0.55 fm. (4) 

Equation (j4]) imposes a numerical lower limit for the pion 
impact parameter, 6, in pionic configurations. Note that 
this number represents a rough estimate, as the inter- 
pretation of RMS radius in terms of a "size" depends 
on the shape of the transverse distribution of partons 
in the core. A more refined estimate, which takes into 
account the intrinsic transverse size of the pion as well 
as the effect of the recoil of the spectator system, is ob- 
tained by requiring that 6/(1 — y) > (i?core + ^V)^^^- 
Assuming that i?J ranges between zero and i?core) ^'^d 
anticipating that the typical y-values in the pion distri- 
bution at 5 ^ i?corc are y = (1 — 2) x Mt^/Mn ~ 0.2, we 
obtain b > 0.44 — 0.62 fm, in good agreement with the 
estimate of Eq. ([4]) . When considering the nucleon's par- 
tonic structure at small x (< 10~^) the above estimate 
of the nucleon core size needs to be modified to account 
for the non-chiral growth due to diffusion in the partonic 
wave function. Also, in this region the transverse size of 
the pion itself can grow due to chiral corrections. These 
effects will be discussed separately in Sees. IVII Al and 
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Chiral dynamics produces also configurations in the 
fast-moving nucleon cfiaracterized by large longitudinal 
separations of the pion and the spectator system, 

I ^ (5) 

with no restriction on b. The relevance of these config- 
urations for the nucleon's partonic structure cannot be 
ascertained without detailed consideration of the effec- 
tive longitudinal sizes of the subsystems and possible co- 
herence effects, and will be discussed in Sec. IVII CI In 
the following we limit ourselves to chiral contributions at 
large transverse distances. 

B. Pion distribution in the nucleon 

In its region of applicability defined by Eqs. (H]) and 
the 6-dependent pion "parton" distribution can be cal- 
culated as the transverse Fourier transform of the "pion 
GPD" in the nucleon. The latter is defined as the transi- 
tion matrix element of the operator measuring the num- 
ber density of pions with longitudinal momentum frac- 
tion y in the fast-moving nucleon, integrated over the 
pion transverse momenta, and with a transverse momen- 
tum transfer Aj^ to the nucleon (see Fig. [2^): 

X (kI 4 J* + A/2) -A/2) 

= {27r)^2P)S^'Hp2~Pi+A)H,{y,t), (6) 

where — P oo, A|| — 0, and 

t ^-Ai. (7) 

Here aj^ and a^^a denote the pion creation and anni- 
hilation operators, and the sum over isospin projections 
(subscript a) is implied. Eq. © refers to the hehcity- 
conserving component of the nucleon transition matrix 
element (A2 = Ai), and H.^{y,t) is the corresponding 
GPD; the helicity-flip GPD is defined in analogously but 
will not be needed in the present investigation. In terms 
of the pion GPD the transverse coordinate distribution 
is then obtained as (& = |6|) 

Uv.h) = J |^e-(^-'')iJ.(y,i). (8) 

We note that a manifestly covariant definition of the pion 
GPD, as the matrix element of a pionic light -ray operator 
between nucleon states, was given in Ref. fli|; the equiv- 
alence of that definition to Eq. ([6]) is shown by going to 
the frame where the nucleon is moving fast and expand- 
ing the pion fields in creation and annihilation operators. 

The pion GPD in the nucleon implies summation over 
all relevant baryonic intermediate states. Because the 
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FIG. 2; The pion GPD in the nucleon. (a) Transition ma- 
trix element of the density of pions with longitudinal momen- 
tum fraction y ~ Mt^/Mm and transverse momentum transfer 
|A^| ~ M^, Eq. ([51 . (b) Invariants used in modeling finite- 
size effects with form factors. ti,2 are the pion virtualities in 
the invariant formulation, Eq. si_2 the invariant masses 

of the TiB systems in the time-ordered formulation, Eq. (|20|l . 

pion wavelength is assumed to be large compared to the 
typical nucleon/baryon radius, only the lowest-mass ex- 
citations can effectively contribute to the GPD in the 
region of Eqs. ((1]) and We therefore retain only the 
N and A intermediate states in the sum: 

Htt = HttN + Ht^a, (9) 

fir = firN + fir A- (10) 

The inclusion of the A, whose mass splitting with the nu- 
cleon introduces a non-chiral scale which is numerically 
comparable to the pion mass, represents a slight depar- 
ture from strict chiral dynamics but is justified by the 
numerical importance of this contribution; cf. the dis- 
cussion of the Nc ^ oo limit in QCD in Sec. IVII 

To study the properties of the 6-dependent pion dis- 
tribution at large distances we need a dynamical model 
which allows us to calculate the pion GPD in the rele- 
vant region of momenta. Here we follow a heuristic ap- 
proach and start from the simplest possible system of 
pointlike pions and nucleons interacting according to a 
phenomenological Lagrangian. We shall see later how 
this definition can be amended to incorporate finite-size 
effects. In which region the results should be regarded as 
physical in the light of the discussion in Sec III AI will be 
the matter of the following investigations. 

The pion GPD in the nucleon can be calculated using 
invariant perturbation theory, by evaluating the matrix 
element in Eq. ([S]), or, equivalently, the matrix element of 
the pionic light-ray operator of Ref. , using the Feyn- 
man rules for pointlike TriV interactions; see Ref. [13] for 
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details. The resulting Feynman integral is computed by 
introducing light-cone coordinates and performing the 
integral over the "minus" (energy) component of the loop 
momentum using Cauchy's theorem. Closing the con- 
tour around the pole of the propagator of the spectator 
baryon, one arrives at a representation in which the spec- 
tator is on mass-shell, and the emitted and absorbed pion 
are off mass-shell, with virtualitics 65] 



ii,2 = fci,2 = -(fc_L T J/Aj^/2)Vj/ + i„ 



(11) 



(see Fig. ^jp) . Here fcj^ is the transverse momentum of 
the spectator baryon. 



y 



(12) 



and 



[yH'll + y{Ml - Ml)] /y (13) 



is the minimum virtuality required by kinematics for a 
given pion momentum fraction, y. The ttN and ttA 
GPDs are then obtained as 

H^N{y,t) = igl^!^ h{y,t-M^,MN), (14) 
if.A(y,t) = 2glj,^Iw{y,t-M^,MA). (15) 

Here gTiNN and QttNA are the coupling constants in the 
conventions of Ref. 14] and Appendix and the dis- 
tributions are the isoscalar pion GPDs, corresponding to 
the sum of 7r+,7r^ and tt'^ distributions in the proton; 
c/. Eq. ([6]). The functions /s and Iiq denote the basic 
transverse momentum integrals arising in the calculation 
of the meson distribution with intermediate octet and 
decuplet baryons, 



_ y f d'^kj^ 



08,10 



Any J {2T:f{ti-Ml){t2-Miy 



(16) 



where 



- [~ti -t2+yt + 2{Mb - MnY] , (17) 



910 



2AMlMl 



[2Ml{-h-t2+t) 



[2 (Ma 



Mi 
+ M 



N) 



h){M 



Ml 



t2 + 



t2)] 

yt] . (18) 



Note that while the ti,2 of Eq. (fTTjl depend on the vector 
A_L, the integral Eq. (fTB|) depends only on t = — 
because of rotational invariance in transverse space. 

As it stands, the transverse momentum integral in 
Eqs. P^ -p5 |) is divergent. This divergence is related to 
short-distance contributions in the pointlike particle ap- 
proximation and does not affect the chiral long-distance 
behavior of the 6-dependent distribution. Several ways 
of regularizing this divergence and extracting the chiral 
contribution will be discussed in the following. 



The pion GPD can equivalently be evaluated in time- 
ordered perturbation theory, where Fig. [5^ is interpreted 
as a process where the fast-moving nucleon (momentum 
P ^ Afjv) makes a transition to a irB intermediate state, 
in which we evaluate the operator measuring the density 
of pions with longitudinal momentum yP, and then back 
to a nucleon state whose transverse momentum differs 
from the original one by Aj_. In this formulation the in- 
termediate particles are on mass-shell, but the energies of 
the ttB states before and after the operator are different 
from that of the initial/final nucleon state. The invari- 
ant masses of the intermediate states, which are directly 
proportional to the energies, are given by 



Sl,2 = (fci,2 +Pb)^ 
^ (fc^T A^/2)^ 

y 



Ml 



kl 



Ml 



y 



(19) 
T- (20) 



(see Fig. [5}d). The connection to the invariant formula- 
tion is established by noting that, for given y and k±, 



A 



Sl,2 



si,2-Mi = 



Ml 



t 



1.2 



y 



(21) 



whence the denominators in Eq. (jl6p can also be inter- 
preted as "energy denominators." The minimum value of 
the invariant mass difference, Asmin, for given momen- 
tum fraction y can be obtained by substituting ti^2 by 
t-min, Eq. (jl3p . Both the invariant and the time-ordered 
formulation will be useful for discussing the properties 
of the chiral long-distance contribution following from 

Eqs. (HSD-KIHI). 



C. Large— 6 asymptotics 

It is instructive to consider the asymptotic behavior of 
the distribution of pions for 6 ^ oo and fixed y. It is 
determined by the leading branch cut singularity of the 
GPD in the i-channel and can be calculated by applying 
the Cutkosky rules to the Feynman graphs of Fig. [2] with 
pointlike vertices |l1J] . The asymptotic behavior is of the 
form 



fiTB{y,b) cx 



-KBb 



KBb 



(22) 



where B = iV, A, . . . denotes the intermediate baryon; 
the expression applies in principle also to higher-mass 
states, cf. the discussion below. The decay constant, kb, 
depends on the pion momentum fraction, y, and is di- 
rectly related to the minimum pion virtuality, Eq. (jl3p , in 
the invariant formulation, or the minimum invariant mass 
difference in the time-ordered formulation, cf. Eq. (|2ip : 



Kb 



(Ml 



y 



1/2 



(23) 



To exhibit the y-dependence of the decay constant in 
the parametric region of chiral dynamics, y M^^/M^^, 
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Eq. 11]), we set 



D. Contribution to nucleon parton densities 



y = tiAU/Mn, 



(24) 



where the scahng variable, 77, is generahy of order unity. 
Substituting Eq. (fT3|) into Eq. ([23]) and dropping terms 
suppressed by powers of M^/Myv, we obtain 



Kb 



[Ml - M%)M^ 



M, 



N 



1/2 



(25) 



This result has several interesting implications: 



(a) For the nucleon intermediate state {B = TV) the 
second term is zero, and one has kn M^r with 
a coefficient of order unity and depending on rj. 
In this case the 6-distribution exhibits a "Yukawa 
tail" with a y-dependent range of the order 

as expected. 

(b) For a higher-mass intermediate state {B ^ N) 
the decay constant is determined by competition 
of the chiral scale, AI^, and the non-chiral scale, 
(M| - MI,)M^/Mn. The larger the N-B mass 
splitting, the smaller 77 has to be for the chiral scale 
to dominate. This effect suppresses the contribu- 
tion of higher-mass baryons to f-K{y,b) at large b 
and finite r]. Note also that the pre-exponential 
factor, which is not shown in Eq. (|^^ for brevity, 
vanishes oc y for y ^ [3] . 

(c) For r] ^ one finds kb 2Mt^ irrespective of the 
N-B mass splitting. In this limit the transverse 
"Yukawa tail" has the range one would naively ex- 
pect from the analogy with the 3-dimensional situ- 
ation. However, this limit is purely formal, as this 
region makes a vanishing contribution to the nu- 
cleon's partonic structure at moderate x; cf. the 
discussion in Sees. Ill Dl and IVlTl below. 

For pion momentum fractions parametrically of order 
unity, y ~ 1, Eq. ([23]) gives a decay constant of the order 
Kb ^ Mjsi. An exponential decay with range ~ l/M^ 
is not a chiral contribution to the pion distribution, as 
is expected, because the values of y lie outside the para- 
metric region of Eq. (jT]). In sum, the large-6 asymptotic 
behavior obtained from the naive pion distribution with 
pointlike ttN couplings fully supports the general argu- 
ments of Sec. Ill Al concerning the parametric region of 
the chiral component. 

One notes that the characteristic transverse range 
of the chiral contribution of the pion distribution, 
l/(2Af7r) — 0.71 fm, is numerically not substantially 
larger than our estimate of the non-chiral "core" size, 
Eq. (ID). This shows that an effective field theory ap- 
proach to chiral dynamics, which implicitly assumes that 
the core has zero size and builds up its structure by 
counter terms, is not practical here, and underscores 
the rationale for our phenomenological approach, where 
finite-size effects are included explicitly. 



The chiral contribution to the nucleon's parton densi- 
ties is obtained as the convolution of the pion momen- 
tum distribution in the nucleon with the relevant par- 
ton distribution in the pion. For the gluon, the isoscalar 
quark/antiquark, and the isovector quark/antiquark den- 
sities it takes the form 1691 



gix, 6)chiral 

dy 



[ — [IttN + frrA] (?/, b) gr,{z), 

Jx y 



[u + d] (x, 6)chiral 

"'dy 

y 



[u + d\ (x, 6)chiral 

[Un + Ua] iy,b) ql°\z), 



(26) 



(27) 



[u - d] (x, 6)chiral 

~'dy 



[d - u] (x, &)chiral 
[lUN-y.A]iy,b)qT\z), (28) 



where 



x/y 



(29) 



is the parton momentum fraction in the pion. Here 
/ttat and /^A are the isoscalar pion distributions (sum 
of 7r+, TT^ and tt^) with N and A intermediate states in 
the conventions of Refs. [l^, [13] and Appendix the 
isovector nature of the asymmetry, Eq. is encoded 

in the numerical prefactors. The functions g^, g*°*, and 
q™' are the gluon, isoscalar (total), and isovector (va- 
lence) quark/antiquark densities in the pion, 

ql°\z) = [u + d\^^^^^{z) = [u + dl^^^^{z) 

= i [u + u + d + d\^^^^^{z), (30) 

q:^\z) = ±[d- u]^^ (z) ^ ±[u^ dl^ (z) 

= ±^[u-u-d + d\^^{zy, (31) 



the latter is normalized as 



dzqT'iz) 



1. 



(32) 



The tt'^ does not have a valence distribution because of 
charge conjugation invariance, and we assume isospin 
symmetry. Note that the parton densities in the pion, as 
well as the result of the convolution integrals in Eqs. (|26|) - 
(PS)) . depend on the resolution scale; we have suppressed 
this dependence for brevity. The convolution formulas 
for the strange antiquark density and the >S'?7(3)-fiavor 
symmetry breaking asymmetry will be given in Sec. IIVI 
The expressions in Eqs. apply to parton mo- 

mentum fractions of the order x ~ Mj^/M^ but oth- 
erwise not exceptionally small, and transverse distances 
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b ^ I/Mtt. In deriving them we have assumed that 
the "decay" of the pion into partons happens locally on 
the transverse distance scale of the chiral 6-distribution, 
b ~ I/Mttt (see Fig. [TJ. This is justified parametrically, 
as for the values of x under consideration the parton mo- 
mentum fraction in the pion does not reach small values 
(x < z < 1 in the convolution integral) and one can ne- 
glect chiral effects which cause the size of the pion itself 
to grow at small z. 

To see in which region of x the chiral contribution to 
the isovector antiquark density is localized, it is conve- 
nient to write the convolution formula Eq. (|28p in the 
form 

a; [J- u] (x, 5) chiral = dy [|/,rAf - j/tta] (y, &) 

J X 

X zq:^\z), (33) 



where we have multiplied both sides of Eq. (|28p by x and 
used Eq. ([29]) on the right-hand side. Now both functions 
in the integrand vanish for small arguments: f-nB^y) ~* 
for y ^ 0, and 2:g™'(z) ^ for z ^ 0. Noting that the 
valence distribution zg™'(z) is localized around z 1/2 
at low scales, and that the pion momentum distribution 
is centered around y ~ Mt^/Mn, we conclude that the 
convolution produces a sea quark distribution in the nu- 
cleon centered around values x — yz ^ (1/2) x Mt^/Mm, 
in agreement with the general expectation. The same 
argument applies to the bulk of the chiral isoscalar den- 
sity, Eq. (|27p , which arises mainly from the valence quark 
content of the pion; only at very small x the non- valence 
quarks in the pion produce a distinct contribution. Note 
also that the valence quark density in the pion at z '--^ 1/2 
is generated mostly by relatively small-size configura- 
tions in the pion, justifying our approximation of ne- 
glecting the intrinsic transverse size of the pion in the 
convolution integrals. 

One immediately sees from Eqs. ([27|) and (I28|l that the 
chiral large-distance component is larger in the isoscalar 
than in the isovector quark distributions, because the N 
and A contributions add in the isoscalar sector, Eq. ((27| . 
while they partly cancel in the isovector sector, Eq. (|28|) 
[l^ . This is contrary to the general expectation that 
chiral effects manifest themselves mostly in the sea quark 
flavor asymmetry d — u. The cancellation between N and 
A contributions in the isovector case becomes perfect in 
the large-iVc limit of QCD and restores the proper Nc 
scaling of the isovector distributions; see Sec. IVll 

In principle one can use the asymptotic expressions 
for the pion distribution in the nucleon, Eqs. ()22p and 
(l23)l . to do a numerical estimate of the large-distance 
contribution to the nucleon parton densities based on 
Eqs. (l26])-(l28l)- This approach was taken in Ref. [14] to 
estimate the chiral contribution to the nucleon's gluonic 
transverse size, proportional to the 6^-weighted in- 

tegral of the impact-parameter dependent gluon density. 
Because of the weighting with this quantity empha- 
sizes large transverse distances, and the estimates of the 



6-integrated chiral contribution are relatively insensitive 
to the lower limit in h imposed in the integral (see also 
Sec. |V|. In the present investigation we are interested 
in the antiquark densities per se (not weighted with 6^), 
where there is no such enhancement of large distances, 
and estimates of the chiral contribution are more sen- 
sitive to the lower limit in b. We therefore approach 
this problem differently, by analyzing the phenomeno- 
logical pion cloud model (which incorporates finite-size 
effects) and establishing down to which b the numerical 
predictions are insensitive to the short-distance cutoff 
(Sec. mil) . The numerical evaluation of the long-distance 
contribution based on Eqs. ((26)) - ([28l) will then be done 
based on the results of this investigation (Sees. HV] and 

El. 



III. PION CLOUD MODEL IN IMPACT 
PARAMETER REPRESENTATION 

A. Modeling finite— size eff'ects 

For a quantitative study of the chiral large-distance 
component in the nucleon's partonic structure we need 
a dynamical model which allows us to compute the dis- 
tribution of pions beyond its leading asymptotic behav- 
ior. In addition, we must address the question down to 
which values of b numerical study of this component is 
meaningful, in the sense that it is not overwhelmed by 
short-distance contributions unrelated to chiral dynam- 
ics. Ultimately, this question can only be answered in a 
dynamical model which smoothly "interpolates" between 
the chiral long-distance regime and the effective short- 
distance dynamics. Here we study this question in the 
framework of the phenomenological pion cloud model, 
where the short-distance dynamics is not treated explic- 
itly, but modeled by form factors implementing a finite 
hadronic size unrelated to chiral dynamics. This study 
serves two purposes — it establishes what part of the 
predictions of the traditional pion cloud model actually 
arises from the long-distance region governed by chiral 
dynamics, and it offers a practical way of computing this 
universal long-distance contribution. 

In the phenomenological pion cloud model, the pion 
GPD in the nucleon is defined by the graph of Fig. [51 cf. 
Eqs. p^ - P^ . in which now form factors are associated 
with the ttNB vertices, rendering the transverse momen- 
tum integral explicitly finite. Two different schemes to 
implement these form factors are commonly used and 
have extensively been discussed in the literature. One, 
based on the invariant formulation in which the spectator 
baryon in on mass-shell, restricts the virtualities of the 
exchanged pions by inserting in Eq. (|16p a form factor 



3^ 



virt 



(34) 



for each nNB vertex (see Fig. [^b) . Here 5'(a) denotes 
a function of finite range which vanishes for a oo; 
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for example, an exponential, exp(— a), or the dipole form 
factor, (1 + a)~^. These form factors can be compared 
to those in the well-known meson exchange parametriza- 
tions of the NN interaction, where the exchanged pion 
is regarded as a virtual particle 9]. The other scheme, 
based on the time-ordered formulation, restricts the in- 
variant mass of the nB systems in the intermediate states 
by form factors of the type 33] 



Sl,2 



Ml 



At 



mv. mass 



(35) 



(see Fig. [^b)- An advantage of this scheme is that 
it preserves the momentum sum rule in the transition 
N —^ ttB, i.e., the longitudinal momentum distribution 
of the baryon B in the nuclcon is given by f-Ksi^ — y) for 
B = N,A m, [111 . The relation between the two differ- 
ent cutoff schemes can easily be derived from Eq. (I^IT) . 
Effectively, 



■'^-virt — 



mv. mass 7 



(36) 



i.e., a constant invariant mass cutoff amounts to a y~ 
dependent virtuality cutoff which tends to zero as y — > 0. 
In the traditional formulation of the pion cloud model, 
without restriction to the largc-6 region, the two schemes 
lead to rather different pion momentum distributions. 
The distributions at large b and y ~ Mt^/Mm, however, 
are dominated by vanishing pion virtualities viz. invariant 
mass differences, so that the results in the two schemes 
become effectively equivalent, up to small finite renor- 
malization effects. In the following numerical studies we 
shall employ the virtuality cutoff as used in Rcf. 10]; the 
equivalence of the two schemes for our purposes will be 
demonstrated explicitly in Sec. IIII CI 

We emphasize that we are interested in the pion cloud 
model with form factors only as a means to identify the 
chiral large-distance contribution and delineate the re- 
gion where it is universal and independent of the form 
factors. We do not consider those aspects of the model 
related to the fitting of data without restriction to large 
distances (tuning of cutoff parameters, -kNB couplings, 
etc.); those have been discussed extensively in the litera- 
ture reviewed in Refs. Id]. 



B. Universality at large h 

We first consider the dependence of the pion distribu- 
tion in the nucleon on the impact parameter, b. Specif- 
ically, we want to demonstrate that it reproduces the 
"universal" chiral behavior Eq. (|22|) at large &, and in- 
vestigate for which values of b the distribution is sub- 
stantially modified by the form factors. To this end we 
calculate the pion GPD by numerical evaluation of the 
loop integral, Eq. ^6]) , with a virtuality cutoff of the type 
of Eq. (j34p . and perform the transformation to the im- 
pact parameter representation according to Eq. ([5]); use- 
ful formulas for the numerical calculation are collected in 



formfactor 
pointlike 

v = 0.3 




b [fm] 

FIG. 3: The transverse spatial distribution of pions in the nu- 
cleon, firNiy, fe), as a function of b, for values y — 0.07 and 0.3. 
Shown is the radial distribution 2Tib fT^Niy, b), whose integral 
over b (area under the curve) gives the pion momentum dis- 
tribution. Solid lines: Pion cloud model with virtuality cutoff 
(exponential form factor, At^n = 1.0 GeV) [l^]. Dashed line: 
Distribution for pointlike particles, regulated by subtraction 
at A.\ = 0; the integral over b does not exist in this case. The 
estimated "core" radius, Eq. ([4]), is marked by an arrow. 



Appendix [B] Figure [3] shows /irAf(y, &) obtained with an 
exponential form factor (Avirt = 1.0 GeV, a typical value 
in traditional applications of the pion cloud model), as 
a function of b for y = 0.07 and 0.3, which is 1/2 and 2 
times Mt^/Mn, respectively. Also shown are the distribu- 
tions obtained with pointlike particles (no form factors), 
in which the loop integral was regularized by subtrac- 
tion at = 0; this subtraction of a A^-independent 
term in the GPD corresponds to a modification of the 
impact parameter distribution by a delta function term 
oc S^'^\b), which is "invisible" at finite b 14]. One sees 
that for 6 > 0.5 fm the results of the two calculations co- 
incide, showing that in this region the pion distribution 
is not sensitive to the form factors. Comparison of dif- 
ferent functional forms of the form factor (exponential, 
dipole) also supports this conclusion. Furthermore, we 
note that for large b both distributions in Fig. [3] exhibit 
the universal asymptotic behavior derived earlier [13] ■ 

It is interesting that the 6-value where in Fig. [3] the 
"universal" behavior of f-!TNiy,b) sets in is numerically 
close to the transverse radius of the nucleon's "core," 
inferred earlier from independent considerations, i?coro ~ 
0.55 fm, cf. Eq. This shows that the pion cloud model 
can safely be used to compute the large-5 parton densities 
over the entire region defined by Eq. 

Figure S] illustrates the connection between the trans- 
verse distance and the pion virtualities in Eq. (|16p from 
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virtuality cutoff, A^j^^ = 1.00 GeV 
inv. mass cutoff, A^j^^ = 1 .66 GeV 



*n = (full) - 




FIG. 4: The median pion virtuality in the unregularized in- 
tegral, Eqs. ([T6|) — (fT8|) . as a function of b, for y — 0.07 (solid 
line) and 0.3 (dotted line). It is defined as the value of the 
virtuality cutoff, A^j^, for which /7rJv(j/,&) reaches half of its 
value for Ajirt — > oo, corresponding to the unregularized in- 
tegral. 



a different perspective. Shown there is the median pion 
virtuahty in the unregularized loop integral, defined as 
the value of the virtuality cutofT, A^j^.^, for which the reg- 
ularized f-^NiU) b) reaches half of its value for A^j^.^. oo; 
the latter coincides with the value obtained by regular- 
ization through subtraction. The function fTiNiy,b) is 
always positive when evaluated with an exponential vir- 
tuality cutoff, and monotonously decreasing as a function 
of A^jj.j , so that the median value of A^j^.^ provides a sen- 
sible measure of the average virtualities in the integral 
Eq. (fT6|l for given y and b. One sees that the average pion 
virtualities in the loop strongly decrease with increasing 
b, indicating the approach to the universal chiral region. 
We recall that the leading asymptotic behavior at 6 — > oo 
is determined by quasi-on-shell pions, cf. the derivation 
in Sec. HTCl 



C. Effective pion momentum distribution 



We now want to investigate the distribution of pions at 
large transverse distances as a function of the momentum 
fraction, y. In keeping with our general line of approach, 
we do this by studying how the momentum distribution 
of the pion cloud model with form factors is modified 
when a restriction on the minimum b is imposed. We 
define the effective momentum distribution of pions with 



FIG. 5: Effective momentum distribution of pions in ttA*' 
configurations with impact parameters b > bo, Eq. (|37p . in 
the pion cloud model. Solid lines: Distributions obtained 
with a virtuality cutoff, Eq. (|34|) (exponential form factor, 
Avirt = 1.0 GeV), for 6o = (full integral), bo = 0.55 fm and 
bo — 1.1 fm. Dashed lines: Same for distributions obtained 
with an invariant mass cutoff, Eq. (|34|l (exponential form fac- 
tor, Avirt = 1.66 GeV). The value of Avirt was chosen such 
that it produces the same total number of pions (y-integral) 
for the full distribution as the given virtuality cutoff. The 
value y = Mtt/Mm is indicated by an arrow. 



6 > &0 as the integral 

j SbQ{b>bo) UB{yM (B = 7V,A); (37) 

for 6o = we recover the momentum distribution of 
pions in the traditional usage of the pion cloud model. 
Figure [S] (solid lines) shows the 6-integrated distribution 
Eq. (j37p . obtained with an exponential virtuality cutoff 
(Avirt = 1.0 GeV), for 6o = (full integral) as well as 
6o ~ 0.55 fm and 1.1 fm, corresponding to 1 and 2 times 
the phenomenological core radius, i?coro- One sees that 
the restriction to large 6-values strongly suppresses large 
pion momentum fractions and shifts the strength of the 
distribution toward values of the order y ~ M^^jM-^, 
in agreement with the general expectations formulated 
in Sec. [Ill From the perspective of the traditional pion 
cloud model, the results of Fig.[5]show that less than half 
of the pions in that model arise from the region b > i?corc, 
where the pion cloud can be regarded as a distinct com- 
ponent of the nucleon wave function 

Also shown in Figure [5] (dashed lines) are the corre- 
sponding distributions obtained with an invariant mass 
cutoff, Eq. . For the sake of comparison the cutoff pa- 
rameter A?jy was chosen here such that it gives the same 
total number of pions (y~integral) for the "full" distri- 
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butions in which no restriction on b is imposed (60 = 0); 
the value of Ajnv — 1.66 GeV thus obtained is within 
the range considered in phenomenological apphcations of 
the pion cloud model [35]. One sees that the full dis- 
tributions are quite different for the virtuality and the 
invariant mass cutoff, as dictated by the relation ([2T|l . 
However, when restricted to large b the y-distributions 
in the two regularization schemes become more and more 
alike, as their strength shifts toward values of the order 
y ~ Mt^/M^. This explicitly demonstrates the equiva- 
lence of the virtuality and the invariant mass regulariza- 
tion in the context of our approach, as announced above. 



D. Extension to 517(3) flavor 

In our studies of the strange sea and the S'C/(3)- 
breaking flavor asymmetry below we shall consider also 
the contributions from K and t] mesons to the sea quark 
distributions at large distances. Because the masses of 
these mesons are numerically comparable to the typical 
hadronic mass scale (as given, say, by the vector meson 
mass) , their contributions to the partonic structure of the 
nucleon cannot be associated with chiral dynamics, even 
at large transverse distances. Still, in the context of the 
present discussion of the pion cloud model, it is instruc- 
tive to study the distribution of K and rj in the impact 
parameter representation, and contrast it with that of 
the TT. 

The pseudoscalar octet meson couplings to the nu- 



0.06 




FIG. 6: Effective momentum distribution of kaons in KA. con- 
figurations with impact parameters b > bo, cf. Eq. (|37|l . in the 
meson cloud model with virtuality cutoflF (exponential form 
factor, Avirt ~ 1.0 GeV). Solid line: bo — (full integral). 
Dashed lines: bo = 0.55 fm and bo = 1.1 fm. 



cleon, as determined by SU{S) flavor symmetry, and 
the definition of their impact parameter-dependent mo- 
mentum distributions are summarized in Appendix [X] 
Significant contributions come only from the KA and 
KJ^* channels. The large-6 behavior of these distribu- 
tions is formally governed by the asymptotic expression, 
Eqs. and with the tt mass replaced by the K 

mass. Figure ini shows the numerically computed effective 
momentum distribution of K in KA configurations, with 
and without restriction to large 6, cf. Eq. ([57)) . which 
should be compared to the corresponding distributions 
for the TT in Fig. [5] One sees that the overall magnitude 
of Jka is substantially smaller than that of /ttN, because 
of the smaller coupling constant (no isospin degeneracy, 
cf. Appendix |X| and the larger meson and intermediate 
baryon mass. More importantly, one notes that the re- 
striction to large b suppresses the K distribution much 
more strongly than the tt distribution; only about 1/5 of 
all kaons in the meson cloud model are located at trans- 
verse distances b > 0.55 fm, and less than 1% are found at 
6 > 1.1 fm. While hardly surprising, these numbers show 
clearly that the K (and 77) contribution to the partonic 
structure above the nucleon's core radius, i?coio — 0.55 fm 
is extremely small. 



IV. LARGE-DISTANCE COMPONENT OF THE 
NUCLEON SEA 

A. Isovector sea d — u 

We now apply the formalism developed in Sees. [Ill and 
mil to study the chiral large-distance contributions to the 
sea quark distributions in the nucleon. To this end, we 
evaluate the convolution formulas, Eqs. ((27l) - ((28| . with 
the 5-integrated pion distribution, Eq. (|37p . where the 
lower limit, bo, is taken sufficiently large to exclude the 
model-dependent small-distance region, cf. Fig. [31 Our 
standard value for this parameter is the phenomenolog- 
ical "core" radius, Eq. (U); variation of this value will 
allow us to estimate the sensitivity of the results to un- 
known short-distance dynamics. While not permitting a 
complete description of the sea quark distributions, our 
results allow us to quantify how much comes from the 
"universal" large-distance region, providing guidance for 
future comprehensive models of the partonic structure. 

We first consider the isovector antiquark distribution 
in the proton, [d — u\{x), which experiences only non- 
singlet QCD evolution and is largely independent of the 
normalization scale. The convolution formula Eq. (|28p 
involves the valence quark distribution of the pion; the 
normalization of this distribution is fixed by Eq. ([^^ , and 
its shape has been determined accurately by fits to the 
ttN Drell-Yan data; see Ref. [S^ and references therein. 
We use the leading-order parametrization of the valence 
distribution provided in Ref. [s^; the differences to the 
next -to-leading order parametrization are minor in this 
case. Figure [3 shows the chiral long-distance contribu- 
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FIG. 7: Solid line: Large-distance contribution to the anti- 
quark flavor asymmetry, [d — u]{x), obtained from tvN and ttA 
configurations restricted to impact parameters b > Rcoic = 
0.55 fm, cf. Eq. ^7}. Dotted lines: Same with b > 0.8 i?corc 
(upper line) and 1.2 i?coro (lower line). Data: Result of anal- 
ysis of final FNAL E866 Drell-Yan data [3]; statistical and 
systematic errors were added in quadrature. All curves and 
data points refer to the scale = 54 GeV^ . 



tion obtained when bo is taken to be the phenomeno- 
logical "core" radius, i?coic — 0.55 fm (solid line), as 
well as the band covered when bo is changed from this 
value by ±20% (dotted lines). Also shown in the fig- 
ure are the results of an analysis of the final data from 
the FNAL E866 Drell-Yan experiment, presented at a 
common scale — 54 GeV^ Q • One sees that the 
large-distance contribution to the asymmetry is prac- 
tically zero for x > 0.3, as expected from the general 
considerations of Sec. [TTl At a; ~ 0.1 the large-distance 
contribution accounts for ^ 30% of the measured asym- 
metry, indicating that most of it results from the nu- 
cleon's core at small transverse distances. This conclu- 
sion is robust and, as demonstrated in Sec. IIIIl does not 
depend on the form factors employed in the calculation 
within the pion cloud model (the specific results shown 
here were obtained with an exponential virtuality cutoff 
with A^N = 1.0 GeV and A^a = 0.8 GeV ^Ifl]). At small 
X 0.01) the large-distance contribution obtained in 
our approach comes closer to the data; however, the qual- 
ity of the present data is rather poor, and it is difficult to 
infer the magnitude of the required "core" contribution 
by comparing the present estimate of the large-distance 
contribution to the data in this region of x. 

One sees from Eq. (^51) that the isovector antiquark dis- 
tribution involves strong cancellations between the con- 
tributions from ttA'^ and ttA intermediate states. This is 
not accidental — the cancellation between the two be- 
comes exact in the large- A^c limit of QCD, and is in fact 



necessary to restore the proper Nc scaling of the isovector 
distribution; see Sec. I VII 

B. Isoscalar sea u + d 

The isoscalar light antiquark distribution, [u + d] (x) , 
is subject to singlet QCD evolution and thus exhibits 
stronger scale dependence than the isovector distribution. 
The convolution formula for this distribution, Eq. (|27p . 
involves the total (singlet) antiquark distribution in the 
pion, which we may write in the form 

ql°\z) = qT\z) + 2ql^^{z), (38) 

where g™' is the valence distribution, Eq. ([3T|) . and q^'^ 
the "sea" distribution [63| 

ql'''' = = = u^- = 4_. (39) 

The pion sea was determined within a radiative parton 
model analysis, supplemented by a constituent quark pic- 
ture which relates the pion to nucleon parton densities, 
and found to be relatively small [s^. Again, we use the 
leading-order parametrization for the parton densities in 
the pion. 

Figure [5] shows our result for the chiral large-distance 
contribution to the isoscalar antiquark distribution, sep- 
arately for the valence and sea distributions in the pion 
as well as the total, at the scale — 2GeV^. One 
sees that the sea in the pion becomes important only at 
x <C Mt^/M}^, where the antiquark momentum fraction 
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FIG. 8: Solid/dashed/dashed-dotted line: Large-distance 
contribution to the isoscalar antiquark density, a;[u-|-d], result- 
ing from -kN and ttA configurations restricted to 6 > /Jcorc = 
0.55 fm. The plot shows separately the contributions aris- 
ing from the valence, sea, and total antiquark density in the 
pion, cf. Eq. Dotted line: MSTW2008LO leading-order 

parametrization [37j . 
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in the pion can become small, 2: <C 1. Altogether, the 
large-distance contribution accounts for only ^ 1/5 of 
the total u + d in the nucleon at a; ~ 0.1. 

The antiquark distribution obtained from nB [B = 
N, A) configurations cannot be larger than the total an- 
tiquark distribution in the nucleon, which includes the 
radiatively generated sea. The large-distance contribu- 
tion calculated in our approach easily satisfies this the- 
oretical constraint, as can be seen from the comparison 
with the parametrization obtained in the recent leading- 
order global fit of Ref. 37] (MSTW20008LO), see FigM 
We note that the traditional pion cloud model without re- 
striction to large 6, which generates pions with transverse 
momenta of the order ~ 1 GeV and virtualities ~ 1 GeV^, 
produces an isoscalar sea which comes close to saturat- 
ing the empirical u + d ai large x with the usual range 
of parameters parameters, and can even overshoot it for 
certain choices [13, [IHl ■ The restriction of ttB configura- 
tions to large b in our approach solves this problem in a 
most natural way. 



C. Strange sea s, s 

The strange sea (s, s) in the nucleon at large distances 
has two distinct components. One is the chiral compo- 
nent, arising from s and s in the pion in ttN and ttA 
configurations. It is given by a similar convolution for- 
mula as the isoscalar sea, u + d, Eq. (|27p . 



0.006 



s(a;, b) 



chiral 



^ dy 



[UN + UA]iy,b)s^{z), (40) 



and similarly for s, where STr{z) and STr{z) are the strange 
(anti-) quark distributions in the pion. Assuming that 
the sea in the pion is mostly generated radiatively [36l |. 
we take them to be equal and proportional to the non- 
strange sea in the pion, Eq. 



(41) 



The other component comes from valence s quarks in 
KY{Y — A, E,I]*) and rjN configurations in the nu- 
cleon. Because the masses of these mesons are numer- 
ically comparable to the typical hadronic mass scale (as 
given, say, by the vector meson mass), their contribution 
to the partonic structure of the nucleon cannot strictly be 
associated with chiral dynamics, even at large transverse 
distances. We include them in our numerical studies be- 
cause (a) it is instructive to contrast their contribution 
to those of ttN and ttA; (b) they contribute to s only 
and could in principle generate different x-distributions 
for s and s, as suggested by the model of Ref. [11] (we 
shall comment on this model below). The couplings of 
the octet mesons to the nucleon, as determined by SU (3) 
symmetry and the quark model value of the F/D ratio, 
as well as the definitions of the corresponding meson mo- 
mentum distributions are summarized in Appendix [XI 
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FIG. 9: Dashed line: Large-distance contribution to the 
strange sea in the nucleon, s = s, from -kN and ttA con- 
figurations (6 > -Rcorc ~ 0.55 fm). Solid line: Large-distance 
contribution to s from KA and K'S* configurations, involving 
the valence strange quark distribution in the kaon; KT, and 
rjN are numerically negligible. Dotted line: MSTW2008LO 
leading-order parametrization of the total s = s [s^], multi- 
plied by 1/10 for easier comparison. 



The contribution of K and rj to s(x, b) in the proton is 
obtained as 



s{x, b) 



— {lfr,N{y,b) Sriiz) 
+ [fKA + .fK^ + .fK^']{y,b)sK{z)}, (42) 



where the factor 2/3 accounts for the probability of the 
J] to be in a configuration with a valence s quark (we 
assume a pure octet state of the 77 and do not take into 
account singlet-octet mixing, as the rj contribution turns 
out to be negligibly small anyway). The functions s^{z) 
and sxiz) are the normalized momentum distributions 
of s in 77 and 



/ dzSrj^K{z) = 1. 
Jo 



(43) 



Assuming SU{^) symmetry, we will approximate these 
distributions by the valence quark distribution in the 
pion. 



Sr).K[z) 



Z). 



(44) 



We again use the leading-order parametrization of 
Ref. [36[ for the valence quark density in the pion. Nu- 
merical evaluation of the meson distributions shows that 
the contributions from rjN and KT, in Eq. ([1^ are negli- 
gible because of their relatively small coupling; we retain 
only the KK and KT* terms in the following. 

Figure [5] shows the different large-distance contribu- 
tions to the strange sea, integrated over b > i?corc = 
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0.55 fm. One sees that for x > 0.1 the large-distance 
strange sea is mostly s coming from the valence s in 
KA and KY,* configurations; the precise magnitude of 
this contribution is sensitive to the lower limit in b, cf. 
Fig. [6l For x < 0.1 the large-distance strange sea in 
the nucleon originates mostly from the strange sea in the 
pion in ttB {B = N, A) configurations, which contributes 
equally to s and s. The different mechanisms result 
in s{x) 7^ s{x) for the large-distance component of the 
strange sea. However, the overall magnitude of the large- 
distance component represents only ~ 1/20 of the empir- 
ically determined average strange sea, ^ [s + s] (x) [37| , 
so that one cannot draw any conclusions about the x- 
distributions of the total s and s in the nucleon from the 
large-distance component. Note that the large-distance 
component at 6 > i?coro represents a much smaller frac- 
tion of the total sea in the case of s and s than for u + d, 
at least in the region x ^ 0.01. 

There are significant differences between the leading- 
order parametrizations of [s-|-s](x) obtained in the global 
fits of Refs. \37] and [39*1; up to a factor '--^ 2 at a; = 0.1. 
However, this does not change our basic conclusion, that 
the large-distance s and s are only a small fraction of the 
total. Also, some of the next-to-leading order fits by sev- 
eral groups ^37„ .40.] have begun to extract information on 
the shapes of s(x) and s{x) individually, by incorporating 
neutrino scattering data which discriminate between the 
two. The difference [s — s] (x) is very poorly determined 
by the existing data, and the fits serve mostly to limit 
the range of allowed values. We note that our approach 
to large-distance contributions and the convolution for- 
mulas of Eqs. remain valid also for next-to- 
leading order parton densities, if the parton densities in 
the pion are taken to be the next -to-leading order ones. 
In the present study we restrict ourselves to the lead- 
ing order, because at this order the parton densities are 
renormalization-scheme-independent and possess a sim- 
ple probabilistic interpretation, and because the present 
comparison of our results with the data does not warrant 
high accuracy. 

We would like to comment on the approach of Ref . [s^ , 
where the shapes of s{x) and s{x) were investigated in 
a light-front wave function model with KA components, 
whose amplitude was adjusted to fit the observed total 
strange sea, [s-|-s](x). As just explained, our results show 
that only a very small fraction of the total s and s sea 
arise from transverse distances b > i?corc ~ 0.55 fm where 
the notion of meson-baryon components in the nucleon 
wave function is physically sensible. Even in the tradi- 
tional meson cloud model without restriction to large &, 
KA configurations with standard form factors [13, IHj 
would account only for ^ 1/4 of the present value of 
s + s (s?! . This shows that the assumption of saturation 
of the strange sea by KA configurations made in Ref. [38| 
would require a KNA coupling ~ 2 times larger than the 
SU{3) value and is not realistic. While we see indica- 
tions for s{x) s{x) in the large-distance contribution, 
and certainly nothing requires the shapes to be equal. 



the magnitude of the effect cannot be reliably predicted 
on the basis of the model of Ref. . 



D. Flavor asymmetry u + d — 2s 

The antiquark SU (3) flavor asymmetry u + d — 2s is 
a non-singlet combination of the isoscalar non-strange 
and strange sea, which exhibits only weak scale depen- 
dence. Since we assume SU{3) flavor symmetry of the 
sea quarks in the pion, Eq. (|41|) . only the valence tt and 
K components of u + d and s enter in this combination 
(we neglect the rjN and KI] contributions): 

[u + d-2s]{x,b) ^ [ ^ [Aiv + ~ 2/ka 

Jx y 

-2fK^,]{y,b) ql^\z). (45) 

The large-distance contribution to this asymmetry is 
shown in Fig. 1101 One sees that the asymmetry over- 
whelmingly results from the valence u and d content of 
the pion in ttN and ttA configuration; the s in the kaon 
of KA and i^E* contributes only at the level of < 10% 
of the pion. Overall, the large-distance contribution ac- 
counts for ~ 1/3 of the observed SU (3) flavor asymmetry 
at X ^ 0.1. 









T ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 

%N + 3tA, , val, h > 




KA + KL* val, Score 




MSTW2008LO 




0^ = 2 GeV^ 


. . . . 1 .... 1 .... 1 .... 1 ... . 



0.1 0.2 0.3 0.4 0.5 



X 

FIG. 10: Solid line: Large-distance contribution to the an- 
tiquark SU{3) flavor asymmetry asymmetry in the nucleon 
u + d — 2s from valence it -I- d in the pion in tvN and nA 
configurations, cf. Fig. [S] (6 > i?coro = 0.55 fm). Dashed 
line: Contribution from the valence s in the kaon in KA and 
A'S* configurations, cf. Fig. |9] Dotted line: Leading-order 
parametrization of Ref. [sj • 
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V. TRANSVERSE SIZE OF NUCLEON 
A. Transverse size and GPDs 

An interesting characteristic of the nucleon's partonic 
structure is the average squared transverse radius of the 
partons with given longitudinal momentum fraction x. It 
is defined as 



{b^)f(x) 



(fbP f{x,b) 



if^q,q,g), (46) 



where f{x,b) is the impact parameter-dependent distri- 
bution of partons, related to the total parton density by 



(fbf{x,b) = fix). 



(47) 



The average is meaningful thanks to the positivity of 
f{x,b) [H, [l3|- Physically, Eq. (|46)) measures the aver- 
age transverse size of configurations in the nucleon wave 
function contributing to the parton density at given x. 
The transverse size implicitly depends also on the scale, 
Q^; this dependence arises from the DGLAP evolution 
of the impact-parameter dependent parton distribution 
and was studied in Ref. 



421 



The average transverse quark/antiquark/gluon size of 
the nucleon is directly related to the i-slope of the cor- 
responding nucleon GPD at t — 0, 



d_ 

di 



Hfix,0) 



(48) 



Here Hf{x,t) = Hf{x,£^ = 0,t) denotes the "diagonal" 
GPD (zero skewness, ^ = 0), with Hf{x, 0) = f{x), which 
is related to the impact parameter-dependent distribu- 
tion as (6 = |b|) 



Hf{x,t 



J dH e-''^^^''^ f{x,b). 



(49) 



For a general review of GPDs and their properties we 
refer to Refs.El. 



B. Transverse size from hard exclusive processes 

By virtue of its connection with the GPDs, the trans- 
verse size of the nucleon is in principle accessible experi- 
mentally, through the t-slope of hard exclusive processes, 

j*{Q^) + N ^ M + N (Af = meson,7, ...), 

at > IGeV^ and \t\ < IGeV^, whose amplitude 
can be calculated using QCD factorization and is propor- 
tional to the nucleon GPDs. In general, such processes 
require a longitudinal momentum transfer to the nucleon 
and probe the "non-diagonal" GPDs (^ ^ 0), so that 
the connection between the observable t-slope and the 



transverse size can be established only with the help of a 
GPD parametrization which relates the distributions at 
^ 7^ to those at ^ = 0. The connection becomes simple 
in the limit of high-energy scattering, ^ « xb/2 <^ 1, 
where the GPDs probed in the hard exclusive process 
can be related to the dia gon al ones in a well-controlled 
approximation; see Refs. j44l . for details. In this ap- 
proximation the amplitudes for light vector meson elec- 
troproduction at small x {(f>, p) and heavy vector meson 
photo/electroproduction (J ftp, T) are proportional to the 
diagonal gluon GPD, and thus 



{da/dt)-''^^^+^ oc H^{x = XB,t). 



(50) 



The gluonic average transverse size can be directly in- 
ferred from the relative i-dependence of the differential 
cross section 



. 4^ 
dt 



da/dt{t) 



d(j/dt{Q) 



1/2 



t=0 



(51) 



The universal t-dependence of exclusive and elec- 
troproduction at sufficiently large and exclusive J /ijj 
photo/electroproduction, implied by Eq. ([50]) . is indeed 
observed experimentally and represents an important 
test of the approach to the hard reaction mechanism; 
see Ref. [i^ for a recent compilation of results. 

Experimental information on the nucleon's gluonic 
size and its dependence on x comes mainly from the 
extensive data on the i-dependence of exclusive J /"ip 
photo/electroproduction, measured in the HERA HI 
[l5| and ZEUS [16|] experiments, as well as the FNAL 
E401/E458 [l7| and other fixed-target experiments; see 
Ref. [431 for a recent summary. The t-dependence of the 
cross section measured in the HERA experiments is well 
described by an exponential, 



{daldt)'^^^-"^+'^ oc exp(Bj/^i), 



(52) 



and assuming that this form is valid near t = 0, the 
nucleon's average gluonic transverse size is obtained as 



(6^ 



2B 



(53) 



For a more accurate estimate, the measured t-slope 
is reduced by '-^ 0.3 GeV^^ to account for the finite 
size of the produced J/tp- The exponential slope mea- 
sured by HI at {W) = 90 GeV is Bj/^ = 4.630 ± 
0.060t;^:?63 GeV"^ [15], and ZEUS quotes a value of 
Bj/^ = 4.15 ± 0.05t^;5^ GeV"^ [11]. The central values 
correspond to a transverse gluonic size at a; ~ 10~^ in the 
range (fe^)g = 0.31 — 0.35 fm'^, substantially smaller than 
the transverse size of the nucleon in soft hadronic interac- 
tions. It is also found that the gluonic size increases with 
log(l/a;) with a coefficient much smaller than the soft- 
interaction Regge slope, c/. the discussion in Sec. IVII Al 
Comparatively little is known about the quark size 
of the nucleon at small x. As explained above, light 
vector meson production at small x couples mainly to 
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the gluon GPD. Interesting new information comes from 
the t~dependence of deeply-virtual Compton scatter- 
ing (DVCS) recently measured at HERA. The HI ex- 
periment 19] obtained an exponential slope of — 
5.45±0.19±0.34 GeV~^ by measuring t through the pho- 
ton transverse momentum; larger by one unit than the 
J/V' slope measured by the same experiment. ZEUS [l^l 
extracted a DVCS slope of = 4.5 ± 1.3 ± 0.4 GeV^ by 
measuring the transverse momentum of the recoiling pro- 
ton, again larger than the J slope measured by that 
experiment; however, the exponential fit to the ZEUS 
data is rather poor and the extracted B^ has large errors. 
We note that in both experiments the B^ values were de- 
termined by an exponential fit over the entire measured 
region of t and thus reflect the average i-dependence, not 
directly the slope at t = 0. Still, the data provide some 
indication that the t-slope of DVCS at t = is larger 
than that of J /ip production (the in the DVCS ex- 
periments here are comparable to the effective scale in 
J/f/' photoproduction, Q^^ w 3GeV^). In leading-order 
(LO) QCD factorization, the DVCS amplitude is propor- 
tional to the singlet quark GPDs, and the t-slope of this 
process is directly related to the nucleon's singlet quark 
size, 

(62),+, = 2B„ (54) 

c/. Eq. ([55)1 . One would thus conclude that 

{b^),+, > (&2)^. (55) 

At next-to-leading order (NLO) the DVCS amplitude 
also involves the gluon GPD, and substantial cancella- 
tion is found between the gluon and singlet quark contri- 
butions to the amplitude. This cancellation amplifies the 
effect of a difference in (6^)^+5 and {b'^)g on the DVCS 
i-slope. Because the gluon contribution is negative and 
cancels ~ 1/2 of the quark contribution [i^, the rela- 
tive change in the slope should be ~ 2 times larger than 
the relative change in the average transverse sizes which 
caused it [49]. 

The quark transverse size of the nucleon in the valence 
quark region is measured in hard exclusive processes at 
Jefferson Lab, in particular with the 12 GeV Upgrade. 
In this kinematics the skewness of the GPDs needs to 
be taken into account (^ ^ 0), and the analysis relies on 
GPD parametrizations. It is interesting that the t-slope 
of p*^ production measured in the recent CLAS experi- 
ment [50i] seems to be compatible with the Regge-based 
GPD parametrization of Ref. [5l| (however, it is presently 
unclear how to describe the absolute cross section within 
this framework). A detailed phenomenological study of 
the transverse distribution of valence quarks, based on 
parton densities and form factor data, was performed in 
Ref. [5J]. 



C. Chiral contribution 

We now want to study the contribution of the chiral 
large-distance region, b ~ 1/M^, to the nucleon's average 
transverse size. Adopting a two-component description, 
we define 

dHb^ [fix, 6)eorc + Q{b > 60) fix, b) 

chiralj 

/?,2\ „ _ J 



= (6^)/, core + (^^^)/, chiral- (56) 

Here fix, 6)coro denotes the parton density arising from 
average configurations in the nucleon, distributed over 
transverse distances b i?corc- The function f{x, 6) chiral 
is the chiral component of the parton distribution, ex- 
tending over distances b ^ 1/Mtt- Following the same 
approach as above, we integrate it over b with a lower 
cutoff, bo, of the order of the core radius, Eq. dH); the 
sensitivity of the results to the precise value of bo will be 
investigated below. Note that in Eq. (f56|) the fe^-weighted 
integral in the numerator is computed in two separate 
pieces, while the denominator in both cases is the total 
parton density (core plus chiral) at the given value of x; 
the (6^) chiral thus defined represents the contribution of 
the chiral component to the overall transverse size of the 
nucleon, not the "intrinsic" size of the chiral component 
alone. 

The "core" contribution to (6^) was estimated in 
Sec. Ill Al and Ref. [l3l, by relating it to the slope of 
the nucleon's axial form factor, which does not receive 
contributions from the pion cloud: 

(6') core ~ 1(^2) axial = 0.3 fm^. (57) 

We have already used this result to fix the short- 
distance cutoff in the integral over the chiral contribu- 
tion. A more quantitative determination of the "non- 
chiral" transverse sizes of the nucleon, including the dif- 
ferences between quarks, antiquarks and gluons and their 
cc-dependence, requires a dynamical model of the nucleon 
which smoothly interpolates between small and large dis- 
tances and will be the subject of a separate study. Here 
we focus on the chiral contribution, (6^)/, chiral, which 
can be calculated in a model-independent manner; we 
compare it to the "generic" core size given by Eq. ([57)) . 
keeping in mind that the latter may have a richer struc- 
ture than reflected by this simple estimate. 

The chiral contribution to the transverse size, Eq. (|56p . 
is obtained by calculating the fe^-weightcd integral of the 
6-dependent pion momentum distribution in the nucleon 
studied in Sec. IIIIl cf. Eq. p7|) . and substituting the re- 
sult in the convolution formula for the nucleon parton 
density, Eq. (j27p et seq. Useful formulas for the numeri- 
cal evaluation of the fe^-weighted integrals are presented 
in Appendix |B] Because of the weighting factor b^ , the 
chiral contribution to the transverse size is much less sen- 
sitive to unknown short-distance dynamics (i.e., to the 
cutoff bo) than the contribution to the parton density 
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FIG. 11: The chiral large-distance contribution to the nu- 
cleon's average transverse size, {b^)f, as defined by Eq. (|56|l . 
as a function of x [Q^ — 3GeV^). Sohd line: Gluonic 
size (/ — g), cf. Ref. p^ . Dashed line: Singlet quark size 
(/ = u + u + d + d) Dotted line: Strange quark size (/ = s + s). 
In all cases, the curves show the sum of contributions from 
ttN and ttA configurations. 



itself, and thus represents a much more interesting quan- 
tity for studying effects of chiral dynamics in the par- 
tonic structure. Furthermore, the 6^-weighted integral 
can reliably be computed using the asymptotic form of 
the distribution of pions at large 6, Eq. ([22|l . as was done 
in Ref. [3] . We can use this to estimate analytically the 
sensitivity of (5^) /,chirai to the lower limit, bo. Evaluating 
the integral 



dHeib>bo) b" Um{vM 



1 



(58) 



with the asymptotic expression Eq. (j22[) . and taking the 
logarithmic derivative with respect to b^, we obtain 

&o dl _ 1 
^1 dh) ^ I 



[y = A'U/Mn), (59) 



where we have used Eq. (|25p for kjv and 6o = ^corc = 
0.55 fm. This shows that the sensitivity of (&^)chirai is 
indeed low — a 20% change in bo causes only a 4% change 

in (5^)/,chiral- 

Our results for the chiral contribution to the nucleon's 
average transverse size and its dependence on x are sum- 
marized in Fig. [Til for the scale = 3 GeV^. The curves 
shown are the sum of contributions from TriV and ttA con- 
figurations; heavier mesons make negligible contributions 
at large distances. For reasons of consistency the nucleon 
parton densities in the denominator of Eq. (j56p were eval- 
uated using the older parametrization of Ref. (ssj , which 
served as input to the analysis of the pion parton dis- 
tributions of Ref. . The following features are worth 
mentioning: 



(a) The chiral contribution to the transverse size is 
practically zero above x ~ Mt^/Mn ~ 0.1 and 
grows rapidly as x drops below this value, in agree- 
ment with the basic picture described in Sec. [Til 
The rise of (6^) /.chiral with decreasing x is more 
pronounced than that of the parton density itself 
because the former quantity emphasizes the contri- 
butions from large distances. 

(b) The singlet u- and d-quark size grows more rapidly 
with decreasing x than the gluonic radius. This has 
a simple explanation: the quark/antiquark density 
in the pion sits at relatively large momentum frac- 
tions z 0.5, while the gluon density in the pion 
requires z < 0.1 to be sizable; because z = x/y 
in the convolution integral, and the pion momen- 
tum fractions are of the order y ~ M^/Mat, the 
relevant values of z are reached much earlier for 
the quark than for the gluon as x decreases below 
M-^/Mn. Thus, the chiral large-distance contribu- 
tion suggests that the transverse quark size of the 



nucleon at 



< 



0.01 is larger than the transverse 
gluon size, cf. Eq. ([551) . The difference between the 
chiral contribution to the average sizes at a; = 0.01 
is 



' q+q, chiral 



'5, chiral 



0.09 fm^ 



(60) 



Assuming identical core sizes for the quark and 
gluon distribution, this would correspond to a dif- 
ference of the leading-order DVCS and J/tp t- 
slopes, cf. Eqs. ^ and ([54l) . 



— B 



1.1 GeV^ 



(61) 



well consistent with the HERA results summarized 
in Sec. IV Bl It should be remembered that the chi- 
ral prediction, Eq. ([CTl) . is for the exact t-slope of 
the cross section a.t t — 0, while the HERA results 
represent effective slopes, obtained by fitting the 
empirical i-dependence over the measured range; 
the comparison may be affected by possible devia- 
tions of the true t-dependence from the exponen- 
tial shape. More quantitative conclusions would 
require detailed modeling of the core contributions 
to the transverse size, which themselves can grow 
with decreasing x due to diffusion, see Sec. I VII XI 

(c) The chiral contribution to the transverse strange 
quark size of the nucleon closely follows that to the 
gluonic size. This is natural, as s -I- s is mostly 
generated radiatively, by conversion of gluons into 
ss, in both the pion and the nucleon. 



VI. PION CLOUD AND LARGE-TV^ QCD 

The relation of the chiral component of the large-6 par- 
ton densities to the large-TVc limit of QCD is a problem 
of both principal and practical significance. First, in the 
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large- A^c limit QCD is expected to become equivalent to 
an effective theory of mesons, in which baryons appear as 
solitonic excitations, establishing a connection with the 
phenomenological notion of meson exchange. Second, our 
calculations show that contributions from A intermediate 
states are numerically large, and the large-iVc limit pro- 
vides a conceptual framework which allows one to treat 
N and A states on the same footing and relate their 
masses and coupling constants. We now want to verify 
that the large-distance component of the nucleon's par- 
tonic structure, calculated from phenomenological pion 
exchange, exhibits the correct A^c^scaling required of par- 
ton densities in QCD (c/. also the discussion in Ref. [ll). 

The general Nc scaling of the unpolarize d q uark den- 
sities in the nuclcon in QCD is of the form [5J] 

gix) - X function(7Vcx), (62) 
[u + d]{x), [u + d]{x) - iV^2 X function (iVcx), (63) 
[u — d]{x), [u — d]{x) ^ iVc X function(A^ca;), (64) 

where the scaling functions are stable in the large-iVc 
limit but can be different between the various distribu- 
tions. Equations (j63p and (j64p were derived by assum- 
ing non-exceptional configurations {x ~ and fixing 
the normalization of the scaling function from the lowest 
moments of the parton densities, i.e., from the condi- 
tions that the total number of quarks scale as Nc, and 
the nucleon isospin as N^. The transverse coordinate- 
dependent parton distributions should generally scale in 
the same manner as the total densities, Eqs. and 
([64)l . as the nucleon radius is stable in the large- A^c limit 
(this applies even to the nucleon's chiral radii, because 

Turning now to the pion cloud contribution to the par- 
ton densities at large b, it follows from the expressions 
of Eqs. (fT4 |) -(fT8 l) and their Fourier transform, Eq. ([H|), 
that the ^-de pen dent distributions of pious in the nu- 
cleon scale as [1J| 

f^N{y,b)Jr,Aiy,b) ^ Nc X iunction{Ncx). (65) 

This behavior applies to values y ~ Mt^/Mn ~ N'^^ 
and values b ~ 7V°, corresponding to \t\ ~ iV° in the 
pion GPD. In arriving at Eq. ([65|) we have used that 

3/2 

Mjv,Ma ~ Nc\ that gT^NN ^ Nc , as implied the 
Goldberger-Treiman relation; and that g-^NA scales in 
the same way as g-^NN- Equation (|65[) states that the 
momentum distribution of pions in the nucleon at large 
Nc scales like that of isoscalar quarks or gluons. At the 
same time, the parton densities in the pion scale as 

gn{z), q-K{z) ^ function(z), (66) 

where z ^ Nc in typical configurations; that is, they 
have no explicit Nc dependence at large Nc. One thus 
concludes that the NcScaling of the convolution integral 
for the pion cloud contribution to the nucleon's antiquark 
densities, for both B ^ N and A intermediate states, is 

/ —UB{y,b)q^{z) ~ N^ xiunction{Ncx). (67) 

Jx y 



This correctly reproduces the general A^c-scaling of the 
isoscalar quark and gluon distribution, Eq. (|63p . where 
the A^ and A contributions are added, cf. Eq. ((28|) . How- 
ever, it may seem that the pion cloud contribution at 
large b cannot reproduce the Nc scaling of the isovector 
distribution, Eq. ()63p . which is suppressed by one power. 
The paradox is resolved when one notes that in the large- 
Nc limit the A^ and A become degenerate, 

Ma-Nn - N-\ (68) 

and their couplings are related by [11] 

g-irNA = ^g-uNN- (69) 

Using these relations one has 

/.A(y,6) = 2Ajv(y,&) {y-N-^), (70) 

as can be seen from Eqs. (fT4 |) - (fT8l) and Eq. keeping 
in mind that t,ti,t2 ~ A^^ the region of interest. By 
virtue of Eq. (|70p the A^ and A contributions at large 
Nc cancel exactly in the isovector convolution integral, 
Eq. ([64|l . ensuring that the result has the proper Nc- 
scaling behavior as Eq. 

In sum, our arguments show that the pion exchange 
contribution at large 6 is a legitimate part of the nu- 
cleon's partonic structure in large-A'c QCD, exhibiting 
the same scaling behavior as the corresponding "aver- 
age" distributions. The inclusion of ttA configurations 
at the same level as ttA^ is essential because they repro- 
duce the proper A^c~scaling of the isovector distributions, 
and because they make numerically sizable contributions 
— twice larger than ttN — to the isoscalar distributions. 

In Ref. [l3| we have shown that the isoscalar large- 
b pion distribution in the nucleon [Jt^n + /tta] {y, b) ob- 
tained from phenomenological soft-pion exchange, can 
equivalently be computed in the chiral soliton picture 
of the nucleon at large Nc, as a certain longitudinal 
Fourier transform of the universal classical pion field 
of the soliton at large transverse distances. Extend- 
ing this connection to the isovector pion distribution, 
[I/ttJV — j/tta] {y,b), which is suppressed in the large- 
Nc limit, remains an interesting problem for further 
study. In particular, this requires establishing the con- 
nection between soft-pion exchange and the collective 
rotations of the classical soliton. 



VII. SMALL a;-REGIME AND LONGITUDINAL 
DISTANCES 

A. Growth of core size through diffusion 

In our studies so far we have focused on chiral contri- 
butions to the nucleon's partonic structure at moderately 
small momentum fractions, x ^ 10~^, which arise from 
individual ttB {B = N, A) configurations in the nucleon 
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wave function. When considering smaller values of x sev- 
eral effects must be taken into account which potentially 
limit the validity of the present approximations. 

One of them is the growth of the transverse size of 
"average" partonic configurations in the nucleon due to 
diffusion. Generally, the partons at small x are decay 
products of partons at larger x] the decay process has the 
character of a random walk in transverse space and leads 
to a logarithmic growth of the transverse area occupied 
by the partons: 

(&^)parton = (6^) parton (^^O ) + '^Ct'^a.rton ln(a;o/a;) 

{x <xq^ 10^2). (71) 

The rate of growth — the effective Regge slope, Oparton 
— depends on the type of parton and generally decreases 
with increasing scale Q^, because higher increase the 
effective transverse momenta in the decay process p^ . 
Measurements of the energy dependence of the t-slope of 
exclusive J/^ production at HERA HI and ZEUS [Ull^ 
indicate that the rate of growth for gluons at a scale 
g2 w SGeV^ is approximately a'g « 0.14 GeV"^ 
significantly smaller than the rate of growth of the trans- 
verse nucleon size in soft hadronic interactions, a'^^^^ ~ 
0.25 GeV^^. Using the former value as a general measure 
of the rate of growth of the nucleon's transverse size due 
to diffusion, we estimate that at ~ 3 GeV^ the trans- 
verse size of the "core" increases from R^^^^ = 0.3 fm^ 
at X = 10-2 to 0.35 (0.4) fm^ at a; = 10-3(10-*). In 
principle this effect pushes the region of ttB configura- 
tions governed by chiral dynamics out to larger b as x 
decreases. However, the rate of growth at this scale is still 
rather small, leaving ample room for such configurations 
in the region x > 10~^. Note that at lower scales the rate 
of growth is larger; studies based on DGLAP evolution 
show that a'g approaches the soft value at ~ 0.4 GeV^ 

El- 

B. Chiral corrections to pion structure 

Another effect which needs to be taken into account 
at small x are modifications of the parton density in 
the pion itself due to chiral dynamics. The same mech- 
anism as discussed above for the nucleon in principle 
works also in the pion itself — the pion can fluctuate 
into configurations containing a "slow" pion and a two- 
pion spectator system. When evaluated in chiral pertur- 
bation theory, the momentum fraction of the slow pion 
relative to its parent in such configurations is of the or- 
der y(7r in tt) ~ Mj^/^AnFj^), where is the pion decay 
constant, and 4:TtFt^ represents the generic short-distance 
scale appearing in the context of the renormalization of 
the loop integrals. Such contributions to the parton den- 
sity and the GPD in the pion were recently computed 
in an all-order resummation of the leading logarithmic 
approximation to chiral perturbation theory (31| . which 
does not require knowledge of the higher-order terms in 



the chiral Lagrangian. For the nucleon parton densities 
this mechanism could become important for x <C 10^^, 
where the effective parton momentum fractions in the 
pion can reach small values z ^ 0.1. In the present 
study we restrict ourselves to nucleon parton densities at 
X ^ 10^2, for which the convolution integrals are dom- 
inated by "non-chiral" values of z. The incorporation 
of such corrections to the partonic structure of the pion 
and extension of the present nucleon structure calcula- 
tion toward smaller x remains an interesting problem for 
future study. In particular, it should be investigated how 
the expressions derived in the leading-log approximation 
of chiral perturbation theory compare to a "single-step" 
calculation of pion structure including finite mass and 
size (form factors), along the lines done here for the nu- 
cleon. 



C. Chiral dynamics at large longitudinal distances 

In our studies in Sees. [TTllVl we considered chiral con- 
tributions to the nucleon's partonic structure at large 
transverse distances, which arise from ttB configurations 
at large transverse separations, b ~ 1/M^. As already in- 
dicated in Sec. Ill Al there is in principle another class of 
ttB configurations governed by chiral dynamics, namely 
those corresponding to large longitudinal separations in 
the nucleon rest frame, 

I ^ 1/M,, (72) 

and arbitrary values of the transverse separation, down 
to 5 = 0. We now want to discuss in which region of x 
such configurations can produce distinct contributions to 
the partonic structure. 

The main limitation in admitting ttB configurations of 
the type Eq. (|7^ as part of the partonic structure arises 
from the possible longitudinal overlap of the relevant par- 
tonic configurations in the pion and the "core." To deter- 
mine the region where this effect plays a role, it is useful 
to consider instead of the parton densities the structure 
fimction for j* N scattering and appeal to the notion of 
the coherence length of the virtual photon. Contributions 
to the partonic structure of the type of the convolution 
integrals of Eqs. correspond to the impulse ap- 

proximation of 7*iV scattering, which requires that the 
coherence length of the process be smaller than the longi- 
tudinal distance between the constituents, so that inter- 
ference effects can be neglected; see e.g. Ref. [56]. Gen- 
erally, the coherence length for 7*iV scattering in the 
nucleon rest frame is given by 

^coh = i2MNx)-\ (73) 

where Mjv is the nucleon mass and x « Q'^/W'^ <C 1 
the Bjorken variable; W is the center-of-mass energy of 
the scattering process. Thus, one would naively think 
that in scattering from a irN system with longitudinal 
separation ~ {2M^)~^ coherence effects set in if a; < 
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Mt^/Mm ~ 0.1. However, this argument neglects the fact 
that in the fast-moving nucleon the pion carries only a 
fraction of the order y ^ M-^/M^ 0.1 of the nucleon's 
momentum, so that the effective center-of-mass energy 
for 7*7r scattering is actually lower by this factor, and the 
coherence length smaller by this factor, than in the 7*iV 
process. Interference effectively takes place only when 
the coherence length for both scattering on the pion and 
on the baryon in the i:B configuration is ^ {2Mt^)~^, 
which requires 

X < 0.01. (74) 

For larger values of x coherence effects are small, and 
there is in principle room for a chiral component of the 
partonic structure at small b and longitudinal distances 
~ 1/M^. In order to calculate this component one would 
need to model the finite-size effects limiting the longi- 
tudinal extension of the pion and the spectator system, 
which is related to the "small-a; behavior" of the par- 
ton densities of the respective systems. We leave this 
problem to a future study. Interestingly, this could re- 
sult in partial "readmission" of the small-6 component of 
the pion cloud model which was excluded in the present 
study, potentially affecting e.g. the comparison with the 
measured flavor asymmetry d — u in Sec. IIV Al 

We note that the interference effects in scattering from 
ttB configurations described here are large in the region 
in which chiral corrections to the structure of the pion 
would become important, c/. the discussion in Sec. lVIIB] 
An interesting question is whether in the chiral perturba- 
tion theory approach these effects come into play alr eady 
at the level of the leading logarithmic approximation [3l| , 
or only at the level of subleading or finite terms. 

VIII. SUMMARY AND OUTLOOK 

The transverse coordinate representation based on 
GPDs represents a most useful framework for study- 
ing the role of chiral dynamics in the nucleon's partonic 
structure. It allows one to identify the parametric re- 
gion of the chiral component {x < AI^t-fM^jb ~ l/M^,-) 
and provides a practical scheme for calculating it in a 
model-independent way. Let us briefly summarize the 
main results of our investigation. 

(a) The contributions from nB {B = N, A) configu- 
rations to the parton distributions become inde- 
pendent of the -kNB form factors at transverse 
distances b > 0.5 fm, and thus can be associated 
with universal chiral dynamics. The lower limit in 
b approximately coincides with the nucleon's core 
radius, i?corc = 0.55 fm, inferred previously from 
other phenomenological considerations. 

(b) Only ~ 1/3 of the measured antiquark flavor asym- 
metry [d — u\{x) Sit X > 10~^ comes from the large- 
distance region b > i?coic, showing that most of 



it resides in the nucleon's core at small transverse 
distances. The traditional pion cloud model, which 
attempts to explain the entire asymmetry from pio- 
nic contributions, gets most of the effect from small 
b where the concept of ttB configurations is not ap- 
phcable. 

(c) The isoscalar antiquark distribution [u + d] (x) ob- 
tained from pions at large b remains safely be- 
low the total antiquark distribution determined by 
QCD fits to deep-inelastic scattering data, leav- 
ing room for the (non-perturbatively and pertur- 
batively generated) antiquarks in the core. This 
naturally solves a problem of the traditional pion 
cloud model, where the pionic contribution can sat- 
urate or even exceed the total antiquark density for 
certain non-exceptional parameter values. 

(d) The strange sea quark distributions, s{x) and s{x), 
overwhelmingly sit at small transverse distances, 
b < i?corc- Neither chiral {ttNjTtA) nor KA con- 
figurations at large b account for more than a few 
percent of the empirical s + s. The predictions of 
Ref. [1^ for the x-dependence of s(a;) and s{x) from 
KA fluctuations rely on the region where the con- 
cept of distinct meson-baryon conflgurations is not 
applicable and require a probability of KA fluctu- 
ations several times larger than what is obtained 
from the standard SU (3) couplings. 

(e) The pionic contributions to the nucleon's transverse 
size, (6^), are much less sensitive to short-distance 
dynamics than those to the parton distributions 
themselves, and thus furnish a new set of clean chi- 
ral observables. The large-distance contributions 
to the nucleon's singlet quark size at x < 0.1 are 
larger than those to the gluonic size, suggesting 
that {b'^)q+q > {b'^)g, in agreement with the pat- 
tern of t-slopes of deeply-virtual Compton scat- 
tering and exclusive J/ip production measured at 
HERA and FNAL. 

In the present study we have limited ourselves to the 
universal large-distance contributions to the partonic 
structure, which are governed by soft pion exchange and 
can be calculated in a model-independent way. A com- 
plete description should include also a model of the short- 
distance part, which actually carries most of the parton 
densities. One way of combining the two would be a 
two-component picture, in which the constituents in the 
"core" act as a source of the chiral pion flelds which prop- 
agate out to distances ~ l/M^r. Such an approach would 
be very effective if the characteristic transverse sizes of 
the "cloud" and the "core" were numerically very differ- 
ent. However, this is not the case — the characteristic 
range of two-pion exchange l/(2M7r) — 0.71 fm is numer- 
ically not much larger than our estimate of the "core" 
size, i?corc = 0.55 fm. Another approach, which appears 
more promising, is based on the idea of a smooth "in- 
terpolation" between the chiral large-distance dynamics 
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and the short-distance regime. In particular, the effec- 
tive theory of Ref. [13], which is based on the large- A^c 
limit of QCD, uses constituent quarks as interpolating 
degrees of freedom; it is valid in a wide region, from dis- 
tances of the order ~ 1/Af^ down to distances of the 
order p k, 0.3 fm — the range of the non-perturbative 
chiral-symmetry breaking forces in the QCD vacuum. It 
leads to a picture of the nucleon as quarks bound by 
a self-consistent pion field (chiral quark-soliton model) 
which is fully field-theoretical and relativistic and 
provides a very good description of the nucleon's quark 
and antiquark densities, including subtle effects such as 
the sea quark flavor asymmetry and its polarization [s^ . 
The results of Rcf. [3] and the present work (see in par- 
ticular Sec. IVip show that this large-iVc description of 
nucleon structure is equivalent to phenomenological soft- 
pion exchange at large transverse distances, thanks to the 
universality of chiral dynamics; it thus, in a sense, con- 
tains the result of the present work as a limiting case. Us- 
ing this large-iVc picture as a script to model the impact 
parameter-dependent parton densities at all h, would cer- 
tainly be an interesting problem for further study. 

Direct experimental study of the chiral component of 
the nucleon's partonic structure through hard exclusive 
processes at x < 0.1 would be possible with a future 
electron-ion collider (EIC). The simplest observables are 
the ^-dependences of the differential cross sections for 
various channels ( J/'0, p, tt) at \t\ <C 0.1 GeV^, and 
their change with x; at sufficiently large such mea- 
surements can be related directly to the t-dependence of 
the gluon and quark GPDs at small t. In particular, such 
measurements should be able to resolve variations of the 
i-slope with t and possible deviations from exponential t- 
dependence. Measurements of exclusive processes require 
high luminosity and the capability to detect the recoiling 
baryon at small angles, which is possible with appropriate 
forward detectors. Another interesting option are pion 
knockout processes, corresponding to exclusive scatter- 
ing from a pion at transverse distances b ^ where 
both the recoiling pion and the nucleon are identified in 
the final state; see Ref. J^] for a detailed discussion. 

The partonic content of the nucleon's pion cloud can in 
principle also be probed in high-energy pp collisions with 



hard processes, such as dijet and Drell-Yan pair produc- 
tion. Such processes, including accompanying spectator 
interactions, are most naturally described in the trans- 
verse coordinate (impact parameter) representation cm- 
ployed in our investigation here. Interesting new effects 
appear in collisions at multi-TeV energies (LHC), where 
the cross sections for hard processes can approach the ge- 
ometric limit (black-disk regime) and the probability for 
multiple hard interactions becomes significant. In this 
situation it is important to realize that the ttB configu- 
rations participate in the high-energy scattering process 
with a fixed transverse orientation, which is frozen dur- 
ing the collision; depending on this orientation one may 
either have a violent collision of the pion with the other 
proton or no interaction at all. The averaging over the 
orientations of the ttB configuration must be performed 
in the colliding pp system with given transverse geome- 
try, not in the partonic wave functions of the individual 
protons. This circumstance affects e.g . the rate of mul- 
tijet events in peripheral collisions [59|. More generally, 
the pion cloud represents an example of transverse cor- 
relations in the nucleon's partonic wave function, which 
are neglected in the usual mean-field approximation for 
high-energy collisions. In particular, such correlations 
play a role in central inclusive diffraction, where they re- 
duce the rapidity gap survival probability relative to the 
mean-field result [6Q] . 
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the transverse center-of-momentum, b, and its distance 
from the center-of-momentum of the spectator system, 
r, can easily be derived for the case of a non-interacting 
system, by starting from the well-known expression for 
the center-of-mass in the rest frame and performing a 
boost to large velocity, taking into account that for the 
non-interacting system the longitudinal momentum frac- 
tions are given by the ratios of the constituent masses 
to the total mass of the system. A more formal deriva- 
tion, based on the light-cone components of the energy- 
momentum tensor, can be found in Ref. p^ . 

[65] In Ref. [ij the pion virtualities were denoted by —s±. 
Here we denote them by ti,2, reserving si,2 for the in- 
variant masses of the ttB system, Eq. (|20p . 

[66] Equation (46) of Ref. flij incorrectly writes the convo- 
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[67] 



[68] 



lution formula for the antiquark flavor asymmetry with 
the antiquark density in the pion rather than the va- 
lence quark density. The correct expression is with the 
valence quark density, cf. Eq. (|28p . This does not affect 
the conclusions about the large-A^c behavior presented 
in Ref. [14], which was the sole point of the discussion 
there. 

The relation of our conventions for the pion par- 



ton densities 
ii;^(GRS) 



to those of Ref. ^ (GRS) is 



The value quoted here corresponds to the arithmetic 
mean of the parametrizations of a'j/^ quoted by the 
HERA HI 15] and ZEUS [3 experiments; see Ref. ^ 
for details. 



APPENDIX A: MESON-BARYON COUPLINGS 
FROM SU{3) SYMMETRY 

In this appendix we summarize the meson-baryon cou- 
plings used in the calculation of SU{3) octet meson 
(tt, K, 7]) contributions to the sea quark distributions (see 
Sec. live]) and derive the expressions for the correspond- 
ing meson momentum distributions in the nucleon. 

For the coupling constants governing the N AI + 
B transitions we rely on SU{3) flavor symmetry, which 
is known to describe the empirical couplings well; see 
Ref. [Ill for a recent review. For 8^8x8 transitions 
there are two independent S'[/(3)-invariant structures, 
with coupling constants traditionally denoted by F and 
D. With the standard assignments of the meson and 
baryon fields the Lagrangian takes the form (we show 
explicitly only the terms describing transitions M + B 
P) 



9s P 



TT n + 



(l-2a)(^i^+I]" 



Tad 



2a)K+K 



3=(1 - 4:a)rip 



V6 



h.c. 



(Al) 



Here = F + D \s the overall octet coupling, which is 
related to our ttNN coupling as 



g-uNN = g-^opp = ffs/v^; (A2) 

we use g-irNN = 13.05 in our numerical calculations p^ . 
The ratio a = F/{F + D) remains a free parameter in 
the context of SU{3) flavor symmetry and can only be 
determined empirically or by invoking dynamical models. 
The SU (6) spin-flavor symmetry of the non-relativistic 
quark model implies F/D = 2/3, and thus 



= 2/5 = 0.4; 



(A3) 



we use this value in our numerical studies in Sec. IIV CI 
For 8 ^ 8x10 transitions there is only a single SU (3)- 
invariant structure, and the Lagrangian is of the form 



^ = 5iop(7j'r+A" 



1 - ^A'°E* 

\/6 V3 



7r-A++ 
h.c. 



(A4) 



The decuplet coupling 510 coincides (up to the sign) with 
our ttAA coupling. 



OttNA 



57r-pA++ — 



-.910; 



(A5) 



we use g-jTNA = 20.22, which is close to the large- 
Nc value of {3/2) g^NN, cf. Eq. Note that our 

definition of the coupling constant g^A^A differs from 
the one of Ref. [3 by a factor, gj^NA (Ref. 0) = 
.97rAfA(this work). 
The GPDs of SU{3) octet meson in the nucleon, 
HMB{y,t)i E^nd the corresponding impact parameter- 
dependent distributions are obtained by straightforward 
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extension of the expressions for ttN and ttA in Sec.|TTl cf. 
Eqs. (0 and (HS]), and Eqs. (ITn])-([THl). We work with the 
isoscalar distributions, in which we sum over the isospin 
components of the intermediate meson-baryon system. 
Using the couphngs provided by Eqs. (|A1[) and (|A4[) . they 
are obtained as (we omit the arguments for brevity) 





— H^ + n + 






4, 




HriN 


— Hfip 




— 3 ^TrNN 


(1- 




Hka 


= Hk+a 




— 3 y-nNN 


(1 + 












(1- 


2a)2 /s, 




















/lO, 








+ Hxos'+ 


- io2 

— 3 i'TrA'A 


/lO. 





(A6) 

Here /§ and /ig are the basic momentum integrals of 
Eqs. p6 l) -(fT8 l) . taken at the appropriate values of the 
meson and baryon masses. Note that because of isospin 
symmetry the distributions for the individual isospin 
components are all proportional to the same function and 
can be expressed in terms of the isoscalar distribution as 



XT 



N, 



etc. (A7) 



where the proportions are determined by the squares of 
the coupling constants in the Lagrangian. The corre- 
sponding 6-dependent distributions, /mb(j/, are then 
obtained by substituting these GPDs in Eq. ([5]). 

To determine the coefficients with which the different 
mesons contribute to a given parton density in the pro- 
ton, one must account for the probability with which the 
parton occurs in the individual meson charge states. For 
example (in abbreviated notation) 

{d- u)p = /^+„ [d- u)^+ + f^op (J - u)^o 

= IUnqT', (A8) 

where 5™' denotes the valence quark distribution in the 
pion, Eq. ([51]) . and we have used Eq. (jA7|) for the im- 
pact parameter distribution and the isospin and charge 
conjugation relations for the parton densities in the pion. 



APPENDIX B: EVALUATION OF 
COORDINATE-SPACE DISTRIBUTIONS 

In this appendix we present expressions suitable 
for numerical evaluation of the transverse coordinate- 
dependent distribution of pions in the nucleon and their 
partial radial integrals. The 6-dependent distribution, 
defined by Eq. ([8]), is evaluated as (we omit the argu- 
ment y and the subscript for brevity) 



fib) 



/ dA^ MA^b) H{t) (B2) 



where b = \b\, A± = |Aj^|, t = — A^, and Jq denotes 
the Bessel function. Equation (|B2p can be used to calcu- 
late the /(&) corresponding to a numerically given H{t), 
as obtained from evaluating the loop integral Eq. (fTB]) 
with ttN form factors. In practice, since H{t) shows only 
a power-like fall-off at large —t > 0, we multiply the 
integrand in Eq. (jB2|) by an exponential convergence fac- 
tor, exp(ei), calculate the integral numerically for finite 
e, and estimate the limiting value for e ^ from the 
numerical data at finite e. 

From Eq. (jB2p we can also derive expressions for the 
partial radial integrals of /(&), including those with a 
weighting factor b^. Using standard identities for inte- 
grals of the Bessel function multiplied by powers of its 
argument, we obtain for the integrals over the region 
b<bo 



dH 6(6 < 60) /(6) - 60 / dA^ Ji(zo) H{t) (B3) 



d'^b e(6 < 60) /(6) 6^ 



dA, 



[2Jo(^o) 



/O ^-L 

+ izo-4/zo)Mzo)]Hit), (B4) 



where 



A_l6o. 



(B5) 



The complementary integrals over the region 6 > 60 are 
calculated by re-writing the original Fourier integral for 
/(6), Eq. jlll), in the form 

j d^b 9(6 > 60) = j dH ~ j d^b 6(6 < 60). (B6) 

The unrestricted integral on the R.H.S. then produces a 
two-dimensional delta function at Aj^ — and can be 
evaluated in terms of H{t ~ 0) or its derivative. In this 
way we obtain 



d% 6(6 > 60) /(6) 



H{t^{)) - I d^be{b<bo) f{b), 



(B7) 



d% 6(6 > 60) 6^ /(6) 
dH 



d^b 6(6 < 60) b' J{b), (B8) 



where the right-hand side can be evaluated using 
Eqs. jm and (|B4|. Note that the A_l integrals rep- 
resenting the 6-integrated distributions, Eqs. (jB3|) and 
(IB4[) . converge more rapidly at large Aj_ than the in- 
tegral representing the original /(6), Eq. (|B2|) . and can 
therefore more easily be computed numerically. 



